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I.  INTRODUCTION 

A.  THESIS  OBJECTIVES 

The  objective  of  this  thesis  is  to  develop  a  method  for  estimating  electron 
densities  and  temperatures  in  the  plasma  jet  of  a  model  magnetoplasma-dynamic  thruster, 
such  as  those  presently  in  development  for  space  propulsion,  by  analyzing  the  broadening 
of  spectral  emissions  within  the  plasma.  Both  Doppler  and  Stark  broadening  effects  are 
analyzed  to  respectively  determine  a  radial  temperature  and  electron  density  profile 
within  the  plasma  arc.  A  computer  program  was  written  to  analyze  the  broadening  of 
hydrogen  line  profiles  by  unfolding  and  examining  the  Stark  and  Doppler  half  widths  of 
the  line  profile.  Three  sets  of  data  taken  at  the  NASA  Jet  Propulsion  Laboratory  (JPL) 
Cathode  Test  Facility  (CTF)  are  analyzed. 

B.  THESIS  OUTLINE 

This  thesis  is  divided  into  six  chapters  and  three  appendices,  describing  the 
background,  conduct,  data  analysis  and  conclusions  of  the  experiment.  Chapter  U 
provides  background  information  for  the  thesis,  giving  a  brief  overview  of  the  basic 
concepts  of  plasmas,  magnetoplasma-dynamic  thrusters,  and  the  theory  of  atomic 
radiation.  Chapter  JJJ  continues  the  background  discussion,  outlining  Doppler 
broadening,  the  Stark  effect  and  Stark  broadening,  and  line  broadening  analysis  of 
spectral  lines.  Chapter  IV  describes  the  experimental  facility  at  NASA  JPL,  as  well  as  the 
data  collection  process.  Chapter  V  describes  the  computer  program  written  to  analyze  the 
data  gathered,  and  presents  and  discusses  the  analysis  of  the  three  data  sets.  Chapter  VI 
presents  conclusions  and  discussion  for  the  furthering  and  improvement  of  this  and  like 
experiments. 

Appendix  A  is  the  listing  of  the  computer  program,  written  in  Interactive  Data 
Language  (IDL),  which  was  written  to  analyze  the  data  collected.  Appendix  B  presents 
the  results  of  the  experiment  and  analysis  in  tabular  form.  Appendix  C  contains  data 
tables  used  for  the  interpolations  in  the  DDL  program. 


II.  BACKGROUND 

A.         ELECTRIC  PROPULSION  DEFINITION 

Electric  propulsion  may  be  defined  as  the  acceleration  of  gases  for  propulsion  by 
electrical  heating  and/  or  by  electric  and  magnetic  body  forces.  As  is  evident  in  the 
definition,  electric  propulsion  may  be  subdivided  into  three  distinct  categories:  (Jahn, 
1968) 

1.  Electrothermal.  Propellant  is  electrically  heated  and  thermodynamically 
expanded  and  accelerated  to  supersonic  speeds  through  a  nozzle,  as  in  a 
chemical  rocket. 

2.  Electrostatic.  Propellant  is  accelerated  by  the  interaction  of  electrostatic 
fields  on  charged  particles,  such  as  ions  or  colloids. 

3.  Electromagnetic.  Acceleration  is  achieved  through  the  interaction  of  electric 
and  magnetic  fields  on  a  highly  ionized  propellant  gas. 

1.         Plasma  Thrusters 

Specifically,  electromagnetic  propulsion  deals  with  the  acceleration  of  a  plasma,  a 
highly  ionized  quasi-neutral  gas,  through  the  interaction  of  currents  driven  through  the 
gas  with  magnetic  fields  established  either  by  those  currents  or  by  external  means  (Jahn, 
1968).  Numerous  ways  exist  of  classifying  electromagnetic  thrusters,  from  thrust  mode 
to  working  fluid  path  geometry,  and  the  names  for  these  various  thrusters  are  equally  as 
varied.  The  work  performed  for  this  thesis  applies  to  a  developmental  facility  for  a 
steady-state,  coaxial,  magnetoplasma-dynamic  (MPD)  thruster. 

a.  The  Coaxial  MPD  Thruster 

A  basic  diagram  of  a  coaxial  MPD  thruster,  showing  its  operation,  is 
shown  in  Figure  2. 1  on  the  following  page.  A  propellant  gas,  commonly  hydrogen,  argon 
or  xenon,  is  fed  into  the  cylindrical  thrust  chamber  and  is  ionized.  Current  is  driven 
through  the  plasma  by  an  external  power  supply  connected  across  the  concentric  anode 
and  cathode.  The  induction  of  a  magnetic  field  by  the  axial  current  component,  or 


the  application  of  a  magnetic  field  externally,  accelerates  the  propellant  by  the  resulting 
(jxB)  force. 

Several  models  of  plasma  thrusters,  both  pulsed  and  steady  state,  have 
been  flown  on  satellites  by  the  United  States  and  the  former  USSR  for  station  keeping 
and  orbit  mainetenance  purposes  (Sutton,  1992).  However,  MPD  thrusters  may  have 
benefits  reaching  far  beyond  their  limited  application  to  date.  Their  high  exhaust 
velocities  and  specific  impulse  make  them  appealing  candidates  for  long  duration,  high 
impulse  space  missions,  such  as  planetary  transfers  and  deep  space  probes. 


Gas  Inlet 


Copper  Anode  (+) 


F,  =  JrB, 


Gas  Inlet 


Figure  2.1.  MPD  Thruster  Operation 

While  some  models  of  plasma  thrusters  have  been  flown,  electromagnetic 
propulsion  is  still  a  largely  developmental  discipline.  MPD  thrusters  employ  both 
electromagnetic  and  plasma  phenomena,  both  of  which  can  prove  very  difficult  to 
analytically  model.  Additionally,  current  designs  are  relatively  inefficient,  and  large 
plasma  thrusters  for  long  duration  missions  have  very  high  power  requirements.  (Sutton, 
1992) 


B.         PLASMA  DEFINITION  AND  PARAMETERS 

Plasma,  the  most  abundant  state  of  matter  in  the  universe,  can  be  simply  described 
as  an  electrified  gas  with  atoms  dissociated  into  positive  (ions)  and  negative  (electrons) 
charged  particles.  The  degree  of  ionization  for  such  a  gas  at  thermal 
equilibrium  is  determined  by  the  Saha  equation: 

^  =  2.4xlO'5I^,-^7  (2.1) 

Here,  n,  and  nn  are,  respectively,  the  number  density  (per  cubic  centimeter)  of  ionized  and 
neutral  atoms,  T  is  the  gas  temperature  in  K,  k  is  Boltzmann's  constant,  and  Ut  is  the 
ionization  energy  of  the  gas  in  question. 

Although  all  gases  have  some  degree  of  ionization,  not  all  gases  may  be  referred 
to  as  plasmas.  A  more  specific  definition  of  a  plasma  is  a  quasineutral  gas  of  charged 
and  neutral  particles  that  exhibits  collective  behavior  (Chen,  1974). 

1.         Plasma  Temperature 

At  thermal  equilibrium,  a  gas  has  a  distribution  of  velocities  and,  therefore, 
kinetic  energies.  Particle  velocities  most  often  obey  a  Maxwellian  or  Gaussian 
distribution.  In  one  dimension,  the  velocity  distribution  is  written  as: 


er    \   j  W        -M2/v,.2    , 

j(u)du  =  — e  du 

vth 


(2.2) 


where 


vth   =  <u 


2kT  (2-3) 


m 


HQT&,f(u)du  is  the  number  of  particles  (per  cubic  centimeter)  with  velocities  between  u 
and  u+du,  and  vth  is  the  thermal  velocity  of  the  particle,  dependent  upon  the  temperature, 
T,  and  mass,  m  of  the  particle. 

The  average  kinetic  energy  of  the  particles  in  the  gas  is: 


2 


*.  -  I  kT  (2.4) 

or,  in  the  three  dimensional  case: 


£  v  =  f  kT  (2.5) 


Equations  2.4  and  2.5  show  the  most  meaningful  representation  of  temperature  as 
it  applies  to  plasmas,  as  a  measure  of  average  energy.  Although  temperature  may  be 
expressed  in  Kelvin,  when  dealing  with  plasmas  it  is  common  to  express  temperature  in 
terms  of  energy,  or  electron  volts  (eV).  Specifically,  a  1  eV  plasma  means  kT  -  1.6xl0"19 
Joules  =  1  eV.  Therefore,  a  1  eV  plasma  has  an  Em  =  1.5  eV  in  three  dimensions. 

2.         Debye  Length  and  Sheaths 

The  requirement  for  quasineutrality  in  a  plasma  implies  that  the  bulk  of  the 
plasma  be  kept  free  of  large  electric  potentials  or  fields.  This  is  accomplished  through  a 
process  known  as  Debye  shielding,  by  which  plasmas  "shield  out"  applied  electric 
potentials. 

The  presence  of  a  large  potential  inside  a  plasma  would  almost  immediately  cause 
a  cloud  of  oppositely  charged  particles  to  surround  the  potential  surface.    This  cloud  of 
charged  particles  is  known  as  a  sheath,  and  it  is  the  presence  of  this  sheath  that  shields  out 
the  large  potential  that  caused  its  formation.  In  a  hypothetical  plasma  with  no  thermal 
motion  (kT  =  0)  the  total  charge  in  the  sheath  would  equal  the  total  charge  in  the 
potential,  and  no  potential  would  exist  outside  the  sheath.  However,  in  a  plasma  of 
temperature  T,  shielding  will  be  imperfect.  At  the  edge  of  the  sheath,  the  energy  of  the 
electric  potential  is  approximately  equal  to  the  thermal  energy  kT  of  the  particles,  which 


could  escape  from  the  potential  well,  allowing  potentials  on  the  order  of  kT/e  to  exist 
within  the  plasma. 

A  measure  of  the  thickness  of  this  sheath,  the  Debye  length,  is  defined  as 


K  = 


\ 


'-  (2.6) 

Anne 


and  is  the  distance  at  which  the  shielded  potential  drops  by  a  factor  of  lie.  Accordingly, 
in  a  collisionless  plasma,  the  field  strength  as  a  function  of  distance  from  the  potential  is 
expressed  as: 


4>  =  (t)ne^D  (2.7) 


Note  that  Te  is  used  in  defining  the  Debye  length,  as  the  electrons  in  the  plasma,  being  in 
general  much  more  mobile  than  the  ions,  do  most  of  the  shielding  by  creating  either  a 
surplus  or  deficit  of  negative  charge. 

The  term  quasineutrality  can  now  be  defined.  In  order  to  be  considered  quasi- 
neutral,  the  effects  of  local  concentrations  of  charge  or  the  introduction  of  external 
potentials  into  the  system  must  be  shielded  out  in  a  distance  short  compared  to  the  overall 
dimensions  of  the  gas.  Doing  so  leaves  the  bulk  of  the  plasma  free  of  large  electric  fields 
and  potentials.  Therefore,  in  order  for  a  gas  to  be  considered  a  plasma,  XD  «  L,  where  L 
is  a  measure  of  the  size  of  the  plasma.  Furthermore,  quasineutrality  implies  that  one  can 
take  ne  *  n,  *  n.  That  is,  ni  and  ne  are  considered  equal,  and  a  common  density,  n,  called 
the  plasma  density,  is  used. 


3.         Plasma  Parameter 

In  order  for  Debye  shielding  to  be  effective  there  must  be  a  sufficient  number  of 
particles  inside  the  sheath.  The  number  of  particles  in  a  "Debye  sphere"  of  radius  AD  can 
be  expressed  as: 


ND  =  ^nnXD  (2.8) 

For  an  ionized  gas  to  be  considered  a  plasma, 

HD  »  1  (2.9) 

4.         Collision  Frequency 

In  order  for  an  ionized  gas  to  behave  as  a  plasma,  the  motion  of  charged  particles 
within  the  plasma  must  be  controlled  by  electromagnetic,  rather  than  hydrodynamic, 
forces.  All  plasmas  have  a  fundamental  plasma  frequency,  expressed  as 


G)     = 

p 


4nne 


\       m  (2.10) 


which  describes  the  oscillation  frequency  of  electrons  inside  the  plasma.  If  tc  is  the  mean 
time  between  collisions  of  charged  particles  with  neutral  atoms,  it  is  required  that 

o)pxc>l  (2.11) 

in  order  for  an  ionized  gas  to  behave  as  a  plasma  rather  than  a  neutral  gas. 
C.         MPD  PLASMA  ARC 

Generating  and  sustaining  the  plasma  arc  in  an  MPD  thruster  entails  the  ionizing 
of  the  propellant  gas  as  well  as  establishing  the  current  and  magnetic  field  necessary  to 
accelerate  the  plasma  jet  to  velocities  high  enough  to  generate  sufficient  thrust.  The 
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propellant  gas,  when  injected  into  the  thruster,  is  immediately  exposed  to  an  ionizing 
potential,  existing  between  the  cathode  and  the  anode.     The  plasma  is  then  accelerated 
past  the  nozzle-shaped  anode. 

1.  Plasma  Arc  Description 

In  steady  state  operation,  the  plasma  is  readily  recognizable  as  an  intensely 
luminous  column  of  gas  extending  from  the  cathode  tip,  called  the  cathode  jet.  Viewed 
externally,  the  exhaust  plume  may  have  a  coaxial  structure,  consisting  of  the  cathode  jet 
and  a  less  intense  coaxial  shell  extending  from  the  lip  of  the  anode  orifice  (anode  jet), 
separated  by  an  even  less  intense  region.    The  whole  plume  usually  extends  many  anode 
diameters  downstream  (Jahn,  1968).  Figure  2.2  shows  a  photograph  of  an  MPD  arc  in 
vacuum  tank  operation.  The  cathode  and  anode  jets  are  easily  distinguishable  exiting  the 
anode  orifice. 


Figure  2.2.  MPD  Thruster  Exhaust  Plume  (Jahn,  1968) 

D.         ATOMIC  STRUCTURE  AND  RADIATION 

The  light  emitted  by  the  plasma  arc  described  above,  and  indeed  in  all  plasmas 
subjected  to  electrical  discharges,  is  limited  to  certain  discrete  wavelengths  that  are 
characteristic  of  the  gases  that  compose  the  plasma.  Light  emission  from  an  atom  occurs 
when  an  electron  within  the  atom  changes,  or  "transitions"  from  an  initial  energy  level  to 


a  lower  energy  level.  The  wavelength  of  the  emitted  photon  is 

X  =  -  (2.12) 

V 

where  c  is  the  speed  of  light  and  v  is  the  frequency  of  the  emitted  photon,  which  is 
directly  related  to  the  change  in  energy  of  the  electron,  as  such: 

E-Ef 

v  =  -i-^  (2.13) 

h 

Here,  Et  and  Ef  are,  respectively,  the  initial  and  final  energy  levels  of  the  electron,  and  h 
is  Planck's  constant.  Furthermore,  numerous  models  formulated  to  describe  atomic 
structure  have  shown  that  the  energy  levels  for  electrons  in  an  atom  are  quantized. 

1.  The  Bohr  Model  for  Hydrogen 

The  simplest  model  of  atomic  structure,  the  Bohr  model,  could  only  accurately 
describe  the  radiation  of  single  electron  atoms  such  as  hydrogen  or  singly  ionized  helium. 
Nonetheless,  the  Bohr  model  accurately  predicts  the  wavelengths  of  all  five  known  series 
of  the  hydrogen  spectrum,  and  works  equally  well  on  singly  ionized  helium.  Although 
certain  postulates  upon  which  the  Bohr  model  is  based  have  since  been  disproven,  the 
model  serves  as  a  suitable  "jumping  off  point"  for  describing  the  hydrogen  spectrum, 
which  is  of  particular  importance  in  this  thesis. 

Bohr  showed  that  the  total  energy  of  an  electron  "orbiting"  the  nucleus  of  an  atom 
of  atomic  number  Z  could  only  take  on  discrete  values 


.  _     mZV      J_  „  .  1A3  (2.14) 


n 


2Tt2     „2 


(4ne0Y2V   n 


Here,  m  is  the  electron  mass,  e  is  the  electron  charge,  and  l/47te0  is  the  Coulomb's  Law 
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constant,  h  is  Planck's  constant  over  2n.  For  hydrogen,  Z  =  1,  and 

En    -    -^T      ^  U    =    1'2'3-  (2.15) 


The  integer  n  is  the  quantum  number,  and  denotes  the  discrete  energy  level  of  the 
electron,  n  =  1,  the  lowest  possible  energy  level,  is  considered  the  ground  state,  while 
higher  energy  levels  are  considered  excited  states. 

2.         Quantum  Numbers 

The  definitions  and  descriptions  put  forth  here  are  by  no  means  complete.  An  in 
depth  treatment  of  quantum  numbers  and  degeneracy  can  be  found  in  Eisberg  and 
Resnick(1985). 

The  time-independent  Schrodinger  equation  relates  the  quantum  state  of  an 
electron  to  its  total  energy: 

-|Iv2i|/(r,e,(|))  +  V(A-)i|;(r,e,(|))  =  ^(r,0,(J))  (2  J6) 

Here,  u  represents  the  reduced  mass  of  the  system,  h  is  Planck's  constant,  and  V(r)  is  the 
Coulomb  potential  of  the  hydrogen  nucleus.  E  is  the  total  energy  of  the  electron,  as 
determined  by  Bohr  for  the  quantum  state  n,  and  is  called  the  eigenvalue. 
The  terms 

-fV2  +  *«  (2.17) 

comprise  an  operator  known  as  the  Hamiltonian,  which  is  designated  H,  such  that 

Hty(r,Q,<b)  =  E^rM)  (2.18) 

i|f(r,6,(j))  is  the  solution  of  the  equation,  called  the  eigenf unction.  It  is  the  product 
of  the  functions  of  the  three  independent  variables  r,  0,  and  (J),  which  describe  the 
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position  of  the  electron  with  respect  to  the  nucleus  in  a  spherical  coordinate  system: 

IW'6^  -  fi„,c-)0,„/e)*„,«t»  (2.19) 

As  is  evident  from  (2. 19),  the  eigenfunctions  are  characterized  by  no  less  than  three 
quantum  numbers,  designated  n,  I,  and  m,.  Possible  values  for  the  three  quantum 
numbers  are  as  follows: 


1'liaZ^-i  (2-20) 

m,  =  -/,-/+!,. ..,0,. ..+/-!,+/ 


The  first  quantum  number,  n,  is  the  quantum  number  from  the  Bohr  model,  and 
describes  the  total  energy  of  the  electron.  En  is  called  the  eigenvalue,  and  n  is  called  the 
principal  quantum  number. 

a.  Orbital  Angular  Momentum 

Each  electron  orbit  associated  with  a  quantum  number  n  has  one  or  more 
associated  values  of  angular  momentum.  The  total  angular  momentum  of  the  system,  L, 
is  quantized  according  to  the  quantum  number  /,  termed  the  azimuthal  quantum  number, 
and  is  expressed  as: 

L  =  y/KlTvjh  (2.21) 


This  only  refers  to  the  total  angular  momentum  of  an  electron  orbiting  a  nucleus.  Lz,  the  z 
component  of  orbital  angular  momentum,  is  quantized  according  to  ml ,  and  may  assume 
the  following  values: 

L>  '  m<*  (2.22) 
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Orbital  angular  momentum  of  the  electron  is  commonly  represented  in  a 
vector  model,  which  displays,  for  a  single  value  of  /,  all  possible  states  of  orbital  angular 
momentum  as  defined  by  m, .  Figure  2.3  below  shows  the  vector  model  of  orbital  angular 
momentum  for  the  state  /  =  2.  Such  vector  models  depict  the  possible  angular  momentum 
states  as  vectors  of  magnitude  L,  which  may  be  located  at  any  point  about  the  z  axis 
defined  by  a  cone  of  height  ±LZ . 


Figure  2.3.  Vector  Model  of  Orbital  Angular  Momentum  for  the  State  /  =  2.  Numbers  are 
in  units  of  K  (From  Eisberg  and  Resnick,  1985) 
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b.  Spin  Angular  Momentum 

In  addition  to  orbital  angular  momentum,  an  electron  posseses  an  intrinsic 
angular  momentum,  S,  called  its  spin.  Spin  is  quantized  according  to  the  quantum 
numbers  s  and  ms,  which  may  assume  the  follwoing  values: 

s  =  1/2 

ms  =  ±1/2  (2.23) 

Spin  angular  momentum  is  quantized  similar  to  orbital  angular  momentum,  specifically: 

s  -  ^T)h 

S    =  m\\ 

z  s 

c.  Total  Angular  Momentum 

The  orbital  and  spin  angular  momenta  of  an  electron  are  coupled  through 
what  is  known  as  the  spin-orbit  interaction.  This  arises  from  an  internal  magnetic  field 
created  by  the  orbiting  of  the  electron  about  the  nucleus.  This  field  is  oriented  in  the 
same  direction  as  L,  and  interacts  with  the  electron's  spin  magnetic  dipole  moment, 
which  is  oriented  in  the  same  direction  as  S,  and  produces  a  torque.  This  torque  serves  to 
couple  the  orbital  and  spin  angular  momenta,  causing  them  both  to  precess  about  their 
vector  sum, 

J  =  L+S  (2.25) 

Here,  J  is  defined  as  the  total  angular  momentum  of  the  electron,  and  is  quantized 
according  to  the  quantum  numbers  j  and  /n  ■ ,  according  to  the  quantization  rules 

Jz  =  mfi  (2.26) 
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where 


J  =  1  +  —  J-  — 

2        2  (2.27) 


m 


\j  =  ~j,  ~j  + 1.. ..,+/-  1,+/ 


Total  angular  momentum  of  the  electron  can  also  be  represented  in  a  vector  model,  much 
like  the  one  used  to  represent  orbital  angular  momentum. 

a.         Degeneracy 

From  (2.20)  it  is  readily  apparent  that  for  every  value  of  n  there  are  several 
different  possible  values  of  /  and  m, .  In  a  one  electron  atom,  such  as  hydrogen,  the  form 
of  the  eigenfunction  depends  on  the  values  of  all  three  quantum  numbers,  but  the  energy 
eigenvalue  depends  only  upon  the  value  of  the  principal  quantum  number,  n.  Because  of 
this,  an  atom  can  have  states  with  completely  different  behavior  as  defined  by  the 
quantum  numbers,  but  have  the  same  total  energy.  This  phenomenon  is  called 
degeneracy,  and  eigenfunctions  corresponding  to  the  same  eigenvalue  are  called 
degenerate.  For  example,  i|/200  and  \}j21±1  describe  different  states  of  a  one  electron  atom, 
yet  each  state  has  the  same  total  energy,  En=2  as  defined  by  the  Bohr  model.  In  reality,  the 
spin-orbit  interaction  results  in  a  minute  difference  in  energy  between  states  with  the 
same  /  but  different./  quantum  numbers.  However,  the  difference  is  very  small,  on  the 
order  of  10 5  eV.  (Eisberg  and  Resnick,  1985) 

Table  2. 1  below  shows,  for  a  one  electron  atom,  all  possible  values  of  the 
quantum  numbers  for  n  -  1 ,  2,  3.  From  Table  2. 1 ,  it  can  be  seen  that  for  every  value  of  / 
there  are  (21  +  1)  degeneracies,  and  for  every  value  of  n  there  are  n2  degeneracies. 
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n 

1 

2 

3 

I 

0 

0 

1 

0 

1 

2 

ml 

0 

0 

-i,o,+i 

0 

-1A+1 

-2,-l,0,+l,+2 

Number  of 
degenerate 
eigenfunctions 
for  each  1 

1 

1 

3 

1 

3 

5 

Number  of 
degenerate 
eigenfunctions 
for  each  n 

1 

4 

9 

Table  2.1.  Possible  Values  of  /  and  m,  for  n  =  1,2,3  (After  Eisberg  and  Resnick,  1985) 

3.         Spectral  Radiation 

An  atomic  electron  will  always  seek  to  occupy  a  state  with  the  lowest  possible 
energy,  within  the  limitations  of  selection  rules  and  the  Pauli  exclusion  principle.  For  the 
hydrogen  atom,  this  energy  level  is  the  ground  state,  defined  by  the  eigenfunction  i|/100 . 
An  atom  in  an  excited  state  will,  over  time,  make  a  series  of  transitions  to  lower  energy 
levels  until  the  ground  state  is  reached. 

Each  transition  to  a  lower  energy  level  is  accompanied  by  the  emission  of  a 
photon  of  energy  £,  -  Ef ,  and  frequency 


E-Ef 


(2.28) 


The  wavelength  associated  with  the  emitted  frequency  is  called  a  spectral  line  or  line,  and 
the  collection  of  lines  resulting  from  all  possible  transitions  within  an  atom  is  called  the 
emission  spectrum  of  the  atom. 

a.  Hydrogen  Spectral  Lines 

The  emission  spectrum  for  hydrogen  can  be  easily  and  accurately 
determined  using  the  Balmer  formula,  which  expresses  the  reciprocal  wavelength  of  each 
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line  as  a  function  of  initial  and  final  energy  states: 


K  =  I  =  Rjyinj       i/n.2) 


(2.29) 


where  /?„  is  the  Rydberg  constant  and  is  defined  as: 

\2.6eV 


R    = 


me 


(47ie0)247t>i3c         2:rhc 


-  1.097xl0'm 


7_  -1 


(2.30) 


b.  The  Balmer  Series 

The  Balmer  series  is  the  family  of  lines  arising  from  transitions  to  nf  =  2 
from  higher  states  (n  =  3,4,5,. ..,n  -  °°).  The  Balmer  series  has  an  infinite  number  of  lines, 
ranging  in  wavelength  from  6562.8  A  in  the  visible  to  the  series  limit  in  -  °°)  3645. 6A  in 
the  near  ultraviolet.  1A  =  1  angstrom  =  10"10  meters.  Table  2.2  shows  some  of  the  more 
prominent  lines  of  the  Balmer  series.  The  Balmer  series  is  one  of  five  known  series  of 
hydrogen  lines,  each  characterized  by  transitions  to  a  particular  nf.  The  other  four  series 
are  the  Lyman  (nf=  1),  Paschen  (nf=  3),  Brackett  (nf=  4)  and  Pfund  (nf=  5)  series. 


Line  Designation 

»« 

A  (A) 

Ha 

3 

6562.8 

HP 

4 

4861.3 

HY 

5 

4340.5 

H6 

6 

4101.7 

He 

7 

3970.1 

Hc 

8 

3889.1 

H„ 

oo  (Series  Limit) 

3645.6 

Table  2.2.  The  Balmer  Series 

The  Balmer  series  is  of  particular  importance  here  because  of  its 
widespread  application  in  the  study  of  Stark  broadening.  Particular  attention  will  be  paid 
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to  the  second  line  in  the  series,  the  Hpline  at  4861.33  A,  as  it  is  the  most  widely  studied 
wavelength  in  the  theory  of  Stark  broadening  (Huddlestone  and  Leonard,  1965). 

c.  Selection  Rules 

Although  the  Balmer  formula  and  the  Bohr  model  accurately  predict  the 
lines  of  all  five  hydrogen  series,  they  are  nonetheless  incomplete  models.  While  they 
accurately  predict  the  change  in  energy  and  wavelength  of  the  resulting  spectral  line,  they 
provide  very  little  insight  into  the  actual  atomic  transition  that  is  taking  place.  Modern 
quantum  theory  has  shown  that  a  transition  can  only  take  place  if  it  satisfies  the  following 
selection  rules: 

A/  -  ±1 

Ay  =  0,±1  (2-31> 

In  order  for  a  transition  to  be  possible  between  two  energy  levels  n,  and  nf ,  their 
respective  /  quantum  numbers  may  only  differ  by  one,  and  their  j  quantum  numbers  can 
only  differ  by  zero  or  one.  Since  the  /  degeneracy  of  the  hydrogen  atom  causes  total 
energy  to  be  only  dependent  upon  the  principal  quantum  number  n,  the  Bohr  model  need 
only  address  changes  in  the  principal  quantum  number  to  successfully  predict  the 
spectrum  of  a  one-electron  atom. 

However,  if  the  /  degeneracy  is  removed,  as  it  is  in  the  presence  of  an 
external  electric  field,  the  total  energy  of  a  quantum  state  becomes  dependent  upon  /,  and 
the  Bohr  model  and  Balmer  formula  can  no  longer  accurately  predict  the  resulting 
spectral  lines.  The  removal  of  the  /  degeneracy  due  to  an  applied  electric  field  is  the 
mechanism  that  causes  the  Stark  effect,  after  which  the  phenomenon  of  Stark  broadening 
takes  its  name. 

Figures  2.4  and  2.5  show  two  different  representations  of  energy  level 
diagrams  for  hydrogen,  and  the  transitions  that  can  result  due  to  the  applicable  theories. 
Figure  2.4  shows  the  energy  levels  and  transitions  that  are  predicted  by  the  Bohr  model, 
as  well  as  the  locations  of  the  corresponding  spectral  lines.  Note  the  nonlinearity  of  the 
wavelength  and  frequency  scale.  Transitions  that  are  predicted  by  the  Bohr  model  but  do 
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not  occur  are  denoted  by  dashed  lines.  Figure  5  shows  some  of  the  allowed  transitions  in 
the  Balmer  (nf=  2)  series,  taking  into  account  selection  rules. 
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Figure  2.4.  Energy  Level  and  Transition  Diagram  for  Hydrogen  According  to  the  Bohr 
Model.  (Eisberg  and  Resnick,  1985) 
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Figure  2.5.  Some  Allowed  Transitions  for  the  Balmer  Series. 
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III.  LINE  BROADENING  IN  PLASMAS 

The  two  major  classes  of  line  broadening  are  Doppler  broadening,  and  collision 
broadening.    While  Doppler  broadening  relates  to  the  velocity,  and  ultimately,  the 
temperature  of  the  radiating  atom,  collision  broadening  is  the  result  of  collisions  on  the 
atomic  level  with  other  particles,  either  charged  or  neutral.  The  two  dominant  causes  of 
line  broadening  in  plasmas  are  Stark  broadening,  due  to  the  perturbation  of  the  energy 
states  of  the  emitting  atom  by  the  electric  fields  of  nearby  charged  particles,  and  Doppler 
broadening,  caused  by  the  thermal  motion  of  the  emitter. 
A.         TERMS  AND  CONVENTIONS 

1.  Line  Shape 

Ideally,  a  spectral  line  is  represented  as  a  line  of  zero  width  at  the  wavelength  at 
which  it  is  emitted.  However,  even  when  completely  isolated  from  its  neighbors  and 
shielded  from  all  external  fields,  a  transitioning  atom  will  emit  a  line  of  finite  width 
(Marr,  1968).  The  natural  width  of  these  lines  tends  to  be  on  the  order  of  hundredths  or 
thousandths  of  an  angstrom,  although  this  is  rarely  the  observed  case.  Typically,  a  line  is 
broadened  by  one  or  more  mechanisms,  and  takes  on  one  of  several  distinct  shapes 
known  as  a  line  shape. 

Line  shape,  or  line  profile,  is  a  term  used  to  describe  the  distribution  of  intensity 
of  a  spectral  line  over  wavelength  or  frequency.  The  shape  is  denoted  /(AA)  or  /(Av), 
where  /  is  the  intensity  distribution. 

AA  =  A-A0 

Av  =  v-v0  (3.1) 

A0  and  v0  are  the  wavelength  and  frequency  of  the  emission,  also  called  the  center 
wavelength  and  frequency.  /(AA)  can  take  on  one  of  three  different  shapes. 
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a.  Lorentzian  Profile 

Also  known  as  a  "dispersion-type"  (Huddlestone  and  Leonard,  1965) 
shape,  the  Lorentzian  shape  is  the  natural  line  shape  of  an  unbroadened  spectral  line 
(Marr,  1968),  and  is  represented  by  the  function 


/(AA) 


1 


1   + 


AX2 


AA 


(3.2) 


1/2 


AA1/2 ,  the  half  width  of  the  profile,  is  the  wavelength  difference  at  which  the  intensity 
drops  by  a  factor  of  two.  2AA1/2  is  known  as  the  full  half-width  of  the  profile.  Lorentzian 
profiles  are  associated  with  collision  broadened  lines,  of  which  Stark  broadening  is  a 
type.  A  lorentzian  profile  is  shown  below  in  Figure  3. 1. 


Figure  3.1.  Lorentzian  Profile,  Axm  =  0.707 
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b.  Gaussian  Profile 

A  Gaussian  shape,  also  known  as  a  Maxwellian  distribution,  takes  the 


form 


-  I0e    AA° 


(3.3) 


where  I0  is  the  intensity  at  the  line  center,  and  XD,  the  Doppler  width,  is  the  wavelength 
difference  at  which  the  intensity  is  reduced  by  a  factor  of  Me.  The  half  width  of  a 
Gaussian  line  shape  is 

(3.4) 
AXm  =  MDJka 

The  term  Doppler  width  is  used  here  in  reference  to  Doppler  broadening,  which  produces 
line  shapes  that  are  distinctly  Gaussian  (Griem,  1964).  Doppler  broadening  will  be 
discussed  in  detail  later  in  this  chapter.  An  example  of  a  Gaussian  profile  is  shown  in 
Figure  3.2. 
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Figure  3.2.  Gaussian  Profile 
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c.  Voigt  Profiles 

In  many  instances,  a  spectral  line  will  exhibit  neither  a  purely  Lorentzian 
nor  a  purely  Gaussian  shape,  but  will  instead  show  features  of  both  distributions.  Such  a 
profile  is  called  a  Voigt  profile,  and  is  formed  by  the  convolution  of  a  Gaussian  with  a 
Lorentzian  profile: 

I(AX)  =  f°°ID(AX)IL(AX*-AX)d(AX)  (3  5) 

The  resulting  profile  is  expressed  as  a  function  of  the  dimensionless  variable  w  (Marr, 
1968): 

/(M)  "  f"    2    '^     ,2  ^  (3-6) 

-°°  nA    +   (w-  v  r 

where 

X-XQ 
W  "  ~Kx;  (3-7) 

and  a  is  called  the  damping  ratio  and  is  defined  as  the  ratio  of  the  Lorentzian  to  the 
Doppler  half  widths: 

AX, 


a  - 


AX^Jha  (3-8) 


*D 


A  Voigt  profile  will  be  the  resultant  line  shape  when  both  collision-type  and  Doppler 
broadening  are  present  in  a  line  profile.  Examples  of  Voigt  profiles  with  varying 
parameters  can  be  seen  in  Figure  3.3. 
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Figure  3.3.  Voigt  Profiles  for  Damping  Parameter  (a)  of  0,  0.5 
B.         DOPPLER  BROADENING 

Doppler  broadening  takes  its  name  from  the  so-called  Doppler  effect,  which  states 
that  radiation  emitted  from  a  source  moving  towards  or  away  from  an  observer  is  shifted 
in  frequency  by 


Av  =  ±v0^ 
c 


(3.9) 


where  v0  is  the  frequency  of  the  source  at  rest,  v,  is  the  velocity  of  the  source,  and  c  is  the 
speed  of  light.  Radiation  from  a  source  moving  toward  an  observer  experiences  a 
positive  shift,  while  a  negative  shift  occurs  when  the  source  is  moving  away  from  the 
observer.  The  resulting  wavelength  shift  is  also  expressed  as 


A*  =  tXq^- 
c 


(3.10) 
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where  X0  is  the  unshifted  wavelength.  In  a  plasma,  the  random  motions  of  radiating 
particles  causes  a  Doppler  shift  in  the  emitted  radiation,  which  manifests  itself  in  a 
phenomenon  known  as  Doppler  broadening.  Assuming  the  motion  of  particles  within  the 
plasma  is  of  a  thermal  nature,  the  velocity  distribution  of  the  emitting  atoms  is  Gaussian, 
and 

dN         1      "—     dvr 


represents  the  fraction  of  particles  moving  in  the  line  of  sight  with  velocities  between  vs 
and  vs  +  dvs ,  and  v  =  (2kT/m)Vl  being  the  most  probable  velocity.  By  substituting  AX  for 
vs  from  (3.10),  and  defining 


a*d  =  ~K  (3.12) 

c 


(3. 1 1)  may  be  rewritten  as: 


^    -  -!—e'^'  d(AX)  (3.13) 


Lastly,  assuming  the  total  intensity,  /,  is  proportional  to  N,  and  I(AX)d(AX)  is  proportional 
to  dN,  the  formula 

I(M)  =  — !_e    Ako  (3.14) 


XDyfii 


is  arrived  at  for  the  line  shape  of  a  Doppler  broadened  line. 
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It  is  apparent  from 


c 


\ 


2kT 

(3.15) 


m 


that  if  the  Doppler  width  of  a  profile  is  known,  the  plasma  temperature,  at  least  for  the 

emitting  atoms,  can  be  determined.  Doppler  broadening  can  also  result  from  the  bulk 

motion  of  the  plasma,  and  is  evident  in  many  pulsed  plasmas  (Huddlestone  and  Leonard, 

1965). 

C.         STARK  BROADENING 

Unlike  Doppler  broadening,  Stark  broadening  is  only  a  single  category  of  what  is 
known  as  collision  or  pressure  broadening.  Collision  broadening  is  largely  due  to  the 
interaction  of  the  emitter  with  surrounding  particles,  and  is  subdivided  into  three  distinct 
types: 

1 .  Resonance  broadening  due  to  interaction  with  atoms  of  the  same  kind. 

2.  Van  der  Waals  broadening  due  to  interaction  with  neutral  molecules  or  atoms 
of  different  species. 

3.  Stark  broadening  due  to  interaction  with  charged  particles  (ions  and  electrons). 

Named  after  the  so-called  Stark  effect,  Stark  broadening  is  the  broadening  of  the  natural 
(Lorentzian)  shape  of  a  spectral  line  due  to  the  interaction  of  the  emitter  with  the 
Coulomb  forces  of  nearby  ions  and  electrons,  which  are  present  in  abundance  in  dense 
plasmas. 

1.         Stark  Effect 

The  Stark  effect  takes  its  name  from  the  investigations  of  Stark  in  1913  on  the 
emission  spectrum  of  hydrogen  in  the  presence  of  an  electric  field  (Marr,  1968).  Stark 
discovered  that  a  splitting  of  spectral  lines  occurs  when  the  radiators  are  in  the  presence 
of  an  electric  field,  due  to  the  removal  of  degeneracies  in  the  unperturbed  eigenvalues  of 
an  atomic  electron.  The  quantum  mechanical  treatment  of  the  Stark  effect  is  rather 
complex  compared  to  other  splitting  mechanisms  such  as  the  Zeeman  effect,  and  is  not 
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needed  in  the  analysis  of  a  spectrum  (Herzberg,  1944).  The  theory  will  therefore  be 
presented  in  its  elementary  form,  with  the  appropriate  references  made  to  quantum 
mechanical  conventions. 

a.  Elementary  Stark  Effect  Theory 

In  the  presence  of  an  electric  field  of  strength  F,  the  magnetic  moment 
associated  with  the  total  angular  momentum,  J,  of  an  atom  is  unaffected.  However,  the 
field  will  polarize  the  atom,  by  displacing  the  positively  charged  nucleus  from  the  center 
of  gravity  of  the  negative  charge  distribution.  This  results  in  the  creation  of  an  electric 
dipole  moment  proportional  to  F,  say  aF,  where  a  depends  upon  the  orientation  of  the 
electron  charge  distribution  to  the  field,  and  is  therefore  a  function  of  J.  This  is  shown  in 
Figure  3.4  where  the  nucleus  N  is  shown  separated  from  the  center  of  negative  charge  S 
(Marr,  1968). 

A  precession  of  J  about  the  electric  field  direction  occurs  due  to  the 
interaction  of  the  electric  dipole  with  the  electric  field,  such  that  the  component  My  in  the 
direction  of  the  field  remains  constant.  The  rate  of  precession  increases  with  increasing 
F,  and  the  energy  change  for  given  angular  momentum  state  J  is  given  by  the  product  of 
the  field  strength  and  the  electric  dipole, 


AE  =  KM  F2 

an       AMyr  (3  16) 


where  AT  is  a  constant  denoting  a  fractional  portion  of  My  (Herzberg,  1944). 

It  can  be  seen  from  (3.16)  that  the  energy  shift  due  to  the  Stark  effect  is 
proportional  to  the  square  of  the  field  strength.  This  is  termed  the  quadratic  Stark  effect. 
Also,  it  is  of  note  that  term  components  that  differ  in  the  sign  of  My  have  the  same  energy, 
as  reversing  the  rotation  of  an  orbiting  electron  (changing  M}  to  -M}  )  does  not  alter  the 
electric  dipole  moment. 
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Figure  3.4.  Vector  Diagram  for  Electric  Field  F  Coupling  for  Atoms.  S  is  the  center  of 
gravity  of  the  electrons  with  resultant  angular  momentum  /  perpendicular  to  the  orbital 
plane.  (After  Marr,  1968) 

b.  Stark  Effect  for  Hydrogen 

Unlike  other  atoms,  hydrogen  exhibits  a  strong  linear  Stark  effect,  that  is, 
the  energy  shift  increases  linearly  with  field  strength.  If  a  perturbing  electric  field  of 
strength  F  is  applied  to  an  atom,  the  potential  of  that  field  is: 


eFz  -  eFrcosQ 


(3.17) 


Time  independent  perturbation  theory  states  that  the  total  energy  of  a  state  n  under  the 
influence  of  a  perturbing  potential  is 


I//'     ,l2 

W    =  E    +  H'       +V  *2 + 

"         "  "•"        "(E,   -  E  f 

v     n  n' 


(3.18) 
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where  H'nn.  is  the  interaction  Hamiltonian  between  energy  states  n  and  n',  defined  as 

H'n,n>    =   fV*v%<W  (3.19) 

(3.19)  defines  the  interaction  between  the  two  states  n  and  n '  due  to  a  perturbing 
potential,  v.  In  most  atoms  subject  to  electric  fields,  the  first  order  term  (H'nn)  of  the 
perturbation  vanishes,  and  the  quadratic  term  must  be  used  to  a  first  approximation. 
However,  if  states  are  degenerate  to  the  extent  that  two  states  with  different  /  quantum 
numbers  have  the  same  energy  (as  in  the  hydrogen  atom)  it  is  possible  to  have  a  nonzero 
first  order  effect  provided  that  A/  =  ±1  and  Am  *  0  (Marr,  1968). 

Both  the  Bohr  theory  and  modern  quantum  mechanics  show  that,  for 
hydrogen,  an  energy  level  n  splits  into  2n  -  1  equidistant  energy  levels,  the  energy 
difference  of  which  increase  linearly  with  field  strength  (Herzberg,  1944). 

2.         Stark  Broadening  Theory 

In  the  presence  of  a  significant  quantity  of  ions  and  electrons,  as  in  a  plasma,  a 
radiator  is  subjected  to  electric  fields  which  consequently  perturb  its  normal  energy  levels 
by  way  of  the  Stark  effect.  Varying  electric  field  strengths  within  the  plasma  cause  a 
spreading  of  the  Stark  levels  in  individual  atoms,  resulting  in  a  statistical  broadening  of 
the  spectral  lines,  known  as  Stark  broadening.  In  dense  plasmas,  Stark  broadening  is  the 
predominant  form  of  line  broadening. 

If  the  perturbing  ion  or  electron  is  treated  as  a  classical  particle,  then  the  duration 
of  each  perturbation  in  a  statistical  distribution  of  perturbing  particles  is  defined  by 

Xs  •=   9/Vav  (3.20) 

where  p  is  the  impact  parameter  of  the  perturbing  particle,  and  v^  is  the  average  velocity 
of  the  perturber.  Due  to  the  large  difference  in  the  velocities  of  the  fast  moving  electrons 
and  slow  ions,  xs  varies  greatly  for  the  particles.  This  resulted,  initially,  in  two  theories  of 
Stark  Broadening,  called  the  impact  approximation  and  the  quasi-static  approximation. 
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a.  Quasi-Static  Approximation 

The  quasi-static  approximation  is  used  in  calculating  the  Stark  broadening 
due  to  ions,  as  xs  is  generally  larger  than  the  average  time  between  collisions  with  ions. 
For  that  reason,  the  motion  of  ions  may  be  neglected  completely,  and  a  static  Stark  effect 
is  assumed,  taking  into  account  the  statistical  distribution  of  electric  fields  at  the  radiating 
atom  (Marr,  1968). 

In  computing  a  line  profile  using  the  quasi-static  approximation,  the 
splitting  of  the  line  resulting  from  the  Stark  effect  is  first  calculated,  which  is  defined  as: 

Av    =  C  FM 

'/     u<r  (3.21) 

where  u  =  1  for  the  linear  Stark  effect,  2  for  quadratic,  and  C(/is  the  Stark  coefficient  for 
transitions  between  initial  state  i  and  final  state /in  a  radiating  atom,  and  is  derived  from 
applying  perturbation  theory  (Griem,  1964).  The  line  splitting  is  then  averaged  over  the 
probability  distribution  of  all  different  electric  microfields  within  the  plasma,  denoted 

*W>  -  jr^J[-f]  =  jrWH(V)  (3-22) 

^0  ^0  ^0 

F0,  called  the  normal  or  Holtsmark  field  strength  (Griem,  1964),  depends  only  upon  ion 
density  and  is  defined  as 


F„  = 


^±±Nm  =  l^xKT9^3    Vlcm 


47I£o"  —  (3.23) 


N  is  the  ion  density  of  the  plasma. 

b.         Impact  Approximation 

The  impact  approximation  for  Stark  broadening  must  take  into  account  the 
infrequent  and  short-lived  interactions  of  the  radiator  with  highly  mobile  electrons,  and 
therefore  cannot  assume  a  time-independent  perturbation  as  in  the  quasi-static 
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approximation.  Determining  an  impact-broadened  line  profile  first  requires  the  solution 
of  a  time-dependent  Schrodinger  equation  of  the  form 

tkML  -  [HA  ♦  vwixw  (3-24) 

at 


where  HA  is  the  Hamiltonian  of  the  unperturbed  radiator,  and  Ve(t)  is  the  potential  of  the 
perturbing  electron.  x(f)  describes  the  state  of  the  emitting  electron,  and  is  defined  as 

X(l)  =  TA(t,0)X(0)  (3.25) 

where  TA(t,0)  is  called  the  time  evolution  operator,  which  relates  wave  functions  at  time  t 
to  those  aU  =  0  (Griem,  1964). 

This  equation  must  be  solved  for  TA  for  all  Ve(t)  produced  by  perturbers  passing 
with  various  velocities,  impact  parameters,  and  times  of  closest  approach.  These 
solutions  are  then  substituted  into  a  correlation  function  which  is  statistically  averaged 
over  all  possible  perturber  trajectories.  Finally,  the  Fourier  transform  of  the  correlation 
function  is  taken  in  order  to  obtain  the  impact  broadened  line  profile.  (Griem,  1964) 

3.         Stark  Profiles 

In  order  to  generate  a  Stark  broadened  line  profile  using  both  the  impact  and 
quasi-static  approximations,  the  quasi-static  theory  is  first  applied,  subjecting  the 
radiating  atom  to  a  static  field  causing  energy  level  shifts  and  line  splitting.  After  that, 
the  electron  broadening  is  superimposed  upon  the  resulting  Stark  patterns  (Huddlestone 
and  Leonard,  1965).  Stark  profiles  have  been  generated  and  tabulated,  most  notably  for 
the  Balmer  series  of  hydrogen,  for  use  in  plasma  spectroscopy.  Stark  broadening 
produces  a  Lorentzian  line  profile. 
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The  tabulated  profiles  themselves  are  one  half  of  a  symmetrical  profile,  and  are 
generally  presented  as  functions  of  S(a)  versus  a.  a  is  known  as  the  reduced  wavelength, 
and  is  expressed  as 

a  =  M  „  _J^i  (3.26) 

F0  2ticF0 


F0  is  the  Holtsmark  field  strength,  dependent  upon  ion  density.  From  S(a),  a  profile  in 
the  wavelength  scale  can  be  obtained  from 

I(X)  =  S(a)^-  =  FQlS(a)  (3.27) 

dk 

Stark  profiles  for  hydrogen  have  been  extensively  calculated  using  numerous 
approximations.  The  most  accurate  Stark  tables  for  hydrogen  have  been  generated  using 
the  "Unified  Classical  Path  Theory",  often  called  the  "Unified  Theory"  of  Stark 
broadening,  developed  by  Vidal,  Cooper,  and  Smith  (1970).  While  adhering  to  the  quasi- 
static  approximation  for  ions,  the  Unified  Theory  of  Stark  broadening  develops  a  more 
complex  model  for  electron  collisions.  This  gives  better  agreement  with  experimental 
data  than  the  impact  theory  does,  particularly  in  the  line  wings  and  transition  region 
between  the  line  wings  and  line  center.  Stark  profiles  for  a  given  spectral  line  are 
normally  presented  in  tabular  format  according  to  electron  density  and  temperature  of  the 
plasma  environment  of  the  radiator.  The  profiles  are  well  known  for  several  lines  of  the 
Balmer  series,  and  have  been  computed  for  other  hydrogen  series  as  well,  most  notably 
the  Lyman  series.  All  computed  profiles  are  normalized  such  that 


/  °S(a)da  =  1  (328) 
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The  most  notable  and  distinguishing  aspect  of  Stark  broadened  line  shapes  occurs 
at  relatively  high  electron  (ion)  densities  (Ne  >  3  x  1014cm"3)  where  the  Stark  effect  due  to 
ions  is  great  enough  that  the  spectral  line  splits.  When  this  happens,  the  line  appears  as 
two  Lorentzian  lines  side  by  side,  which  overlap  at  the  line  center.  Figure  3.5  shows  the 
comparison  of  theoretical  and  experimental  Stark  profiles  for  the  Hp  line  at  an  electron 
density  of  6.4xl016  cm"3.  Generally,  the  splitting  of  the  line  is  stronger  in  theory  than  in 
experiment  (Vidal,  Cooper  and  Smith,  1972),  although  in  Figure  3.5  it  is  quite  evident  in 
both  the  measured  and  theoretical  profiles. 
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Figure  3.5.  Comparison  of  Experimental  and  Theoretical  Stark  Profiles  of  the  Hp  Line. 
Ne  =  6.4xl016  cm'3,  Te  =  12,210  K.  (From  Huddlestone  and  Leonard,  1965) 

4.  Electron  Densities 

The  most  common  use  of  Stark  broadening  as  a  diagnostic  tool  is  in  determining 
electron  densities  in  plasmas  (Huddlestone  and  Leonard,  1965).  The  most  common  way 
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of  doing  this  is  by  comparing  the  full  half  width  of  a  spectral  line  with  the  half  widths  of 
computed  Stark  profiles.  Electron  density  is  then  determined  using  the  relation 

A\m  =  2F0am  =  2.5xlO-9oc1/2<3       A  (3  29) 

where  AA1/2  is  the  full  half  width  of  the  line  profile,  and  am  is  the  value  of  a  at  which 
S(a)  is  one  half  its  maximum  value.  am  must  be  retrieved  from  the  theoretical  profile. 
This  has  already  been  done  for  various  hydrogen  lines,  at  various  electron  densities  and 
temperatures,  by  Griem  (1974).    Ne  is  the  electron  density  of  the  plasma.  Using  half 
width  comparison,  errors  on  the  order  of  10  -  15%  can  be  expected,  the  lower  limit  being 
achievable  for  lower  (  1015  cm"3)  electron  densities  (Griem,  1964). 

The  second  way  of  determining  electron  densities  is  by  a  least  squares  fit  of  an 
observed  line  profile  to  the  theoretical  Stark  profiles.  This  method  is  the  better  of  the  two 
for  higher  electron  densities  (~  1016  cm"3)  where  splitting  of  spectral  lines  (especially  the 
Hp  line)  becomes  significant.  At  lower  electron  densities  (  1015  cm"3),  error  should  be 
on  the  order  of  5%,  while  an  upper  limit  of  10%  error  can  be  expected  for  higher  electron 
densities  (Vidal,  Cooper  and  Smith,  1972). 
D.         ELECTRON  DENSITY  DETERMINATION  IN  PLASMAS 

Very  rarely  will  a  spectral  line  emitted  in  a  plasma  be  either  purely  Doppler-  or 
Stark  broadened.  In  fact,  this  is  almost  never  the  case;  both  forms  of  line  broadening  will 
invariably  contribute  to  the  line  shape,  and  a  Voigt  profile  will  be  the  result.  Before  Stark 
broadening  analysis  can  be  applied  to  the  spectral  line,  the  Voigt  profile  must  be  analyzed 
to  determine  the  extent  to  which  each  broadening  mechanism  is  present,  and  the  Doppler 
broadening  must  be  factored  out  of  the  line  profile. 

1.         Half  Width  Estimation 

When  concerned  only  with  the  half- width  of  the  observed  line  profile,  it  is 
sufficient  to  determine  the  Voigt  parameters  which  best  fit  the  observed  profile,  from 
which  the  Lorentzian  and  Gaussian  half-widths  can  be  extracted,  and  density  and 
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temperature  can  be  determined.  The  relationships  between  the  Stark,  Doppler,  and 
observed  half-widths  of  line  profiles  are  presented  in  Appendix  C. 

Once  the  Stark  half  width  of  the  profile  is  known,  it  can  be  compared  to  the 
known  half  widths  of  theoretical  Stark  profiles,  which  are  tabulated  as  functions  of  Ne  and 
T  in  Appendix  HI  of  Griem  (1974)  for  several  hydrogen  and  ionized  helium  lines.  Given 
a  reasonable  value  of  T,  either  measured  or  determined  from  the  Doppler  width  of  the 
profile,  and  the  Stark  width  of  the  observed  profile,  an  electron  density  can  be 
interpolated  from  the  half  width  tables  presented  in  Griem  (1974).  Figure  3.6  shows,  for 
Hp,  electron  density  versus  half  width  for  7=  5000  and  40000  K  . 
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Figure  3.6.  Electron  Density  Versus  Full  Half  Width  of  the  H«  Line  for  T=  5000  and 
40000  K. 
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IV.  EXPERIMENT  SETUP 

A.         CATHODE  TEST  FACILITY 

The  cathode  test  facility  (CTF)  at  NASA-JPL  is  designed  to  observe  and  analyze 
cathodes  under  the  simulated  operating  conditions  of  an  MPD  thruster.  It  contains  all  the 
necessary  hardware  and  software  to  observe  and  record  line  profiles  emitted  from  the 
plasmas  generated  in  the  facility.  A  complete  description  of  the  CTF  is  presented  in 
Goodfellow(1996). 

1.         Vacuum  Chamber 

The  main  component  of  the  cathode  test  facility  at  NASA-JPL  is  a  stainless  steel 
vacuum  chamber,  0.5  m  in  diameter  and  2.4  m  long,  consisting  of  four  water-cooled, 
cylindrical  segments.  A  water-cooled,  copper  liner  is  inserted  in  the  middle  two  segments 
to  facilitate  long  duration  cathode  operation. 

The  cathode  fixture  is  mounted  on  the  vacuum  chamber  door,  and  consists  of  two 
coaxial  tubes,  electrically  isolated  from  each  other  and  the  door.  The  cathode  itself  is 
clamped  to  the  inner  tube,  which  serves  as  the  current  feed  for  cathode  operation.  The 
cathode  protrudes  through  a  water  cooled  copper  disk  mounted  on  the  end  of  the  outer 
tube.  The  propellant  gas  is  injected  between  the  two  tubes,  and  enters  the  discharge 
chamber  (the  first  segment  of  the  vacuum  chamber)  through  an  annulus  formed  around 
the  base  of  the  cathode.  The  cathode  itself  is  76  mm  in  length  and  9.5  mm  in  diameter 
with  a  hemispherical  tip,  and  is  made  of  2  percent  thoriated  tungsten. 

The  anode  is  a  water  cooled  copper  ring  7.6  cm  in  diameter,  which  is  mounted  on 
a  flange  between  the  first  two  segments  of  the  vacuum  chamber,  from  which  it  is 
electrically  isolated.  The  last  segment  of  the  vacuum  chamber  contains  a  heat  exchanger 
composed  of  water-cooled,  finned  copper  tubing,  which  cools  the  exhaust  before  it  enters 
the  pumping  system.    Several  optical  ports  line  both  sides  of  the  chamber,  providing 
access  to  the  discharge  chamber  (particularly  the  cathode  region  itself)  and  the  plasma 
plume.  The  rear  of  the  chamber  has  a  port  which  allows  viewing  of  the  cathode  tip  and 
plume  along  the  axis  of  the  tank. 
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The  chamber  is  pumped  down  using  a  610  liter/sec  Roots  blower,  backed  by  a 
140  liter/sec  Stokes  mechanical  pump.  Vacuums  of  less  than  0.13  Pa  (0.975  mTorr)  can 
be  achieved  without  propellant  flow  into  the  chamber,  while  a  vacuum  of  approximately 
80  Pa  (0.6  Torr)  can  be  achieved  with  an  argon  flow  rate  of  0.75  g/s  (Goodfellow,  1996). 
Typical  operating  pressures  during  the  experiment  ranged  from  1 1  Torr  ( 1 .467  kPa)  to  30 
Torr  (4  kPa).  Pressure  in  the  discharge  chamber  was  controlled  either  by  adjusting 
propellant  gas  flow  rate  into  the  chamber,  or  by  adjusting  a  gate  valve  between  the 
chamber  and  the  pumps. 

Power  to  the  arc  is  supplied  from  two  Miller  welding  power  supplies,  each 
capable  of  supplying  a  continuous  1 500  A  at  60  V,  or  2000  A  at  50  V  for  up  to  20 
minutes.  The  initial  arc  breakdown  is  achieved  with  a  4  A,  850  V  start  supply. 

2.  Tank  Instrumentation 

Arc  current  is  monitored  with  precision  shunts  replacing  the  factory  shunts  on  the 
Miller  power  supplies.  Terminal  voltage  is  measured  at  the  current  feedthroughs  into  the 
vacuum  tank. 

A  Sierra  Instruments  Side-Trak  Model  830  flow  meter  and  a  Micromotion  Model 
D6  flow  meter  are  used  to  monitor  propellant  flow  into  the  thrust  chamber.  Propellant 
flow  is  controlled  with  a  throttling  gate  valve  located  upstream  of  the  inlet  to  the  cathode 
fixture.  An  MKS  Baratron  capacitance  manometer  (0  -  133  Pa  range)  is  used  to  monitor 
tank  pressure. 

The  above  parameters,  as  well  as  various  system  temperatures  (cathode,  anode, 
etc.)  were  monitored  and  recorded  with  a  Macintosh  computer  utilizing  Lab  View 
software  and  Opto-22  data  acquisition  hardware. 

3.  Data  Collection  Equipment 

A  system  for  performing  emission  spectroscopy  within  the  plasma  plume  was 
already  in  place  in  the  cathode  test  facility.  A  lens  was  mounted  outside  an  optical  port  at 
the  thrust  chamber,  which  focussed  an  image  of  the  cathode  on  a  screen.  The  inlet  to  a 
length  of  100  micron  diameter  fiber  optic  cable  was  located  at  the  center  of  the  screen, 
flush  with  the  surface.  The  image  screen  and  fiber  inlet  were  both  mounted  on  a 
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micrometer-operated  translation  stage,  which  allowed  the  fiber  inlet  to  be  placed  at  any 
point  within  the  image.  The  translation  stage  allowed  adjustments  in  position  to  less  than 
0.5  mm. 

Figure  4.1  shows  a  basic  schematic  of  the  data  collection  system.  Light  exiting 
the  fiber  was  focused  on  the  entrance  slit  of  a  one  meter  McPherson  scanning 
monochromator.  A  1200  groove/mm  grating  was  used  in  the  monochromator  to  disperse 
the  light,  and  the  variable  entrance  and  exit  slits  were  typically  set  at  30  microns. 


Image  Plane 


Vacuum  Chamber 


Optic 


Monochromator 


Amplification 


Figure  4. 1 .  CTF  Data  Collection  Equipment  Setup. 

Spectral  resolution  of  the  spectrometer  is  on  the  order  of  1  A.  The  output  of  the 
monochromator  is  amplified  and  stored  as  a  voltage  in  a  Macintosh  computer,  again  using 
LabView  software. 
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B.         DATA  COLLECTION 

Data  were  gathered  while  the  CTF  was  operating  on  an  argon  plasma,  with  small 
amounts  of  hydrogen  (10  to  15%  by  volume)  intermixed  to  allow  observation  of  the  Hp 
line.  Under  normal  operating  conditions,  the  plasma  arc  flows  from  the  tip  of  the 
cathode,  forming  a  cathode  jet  as  defined  in  Chapter  n.  Sets  of  hydrogen  line  profiles  for 
analysis  were  gathered  from  within  this  cathode  jet.  The  scanned  wavelength  interval 

o  o 

was  from  4856A  to  4866A,  to  allow  inclusion  of  as  much  of  the  wings  of  the  line  as 
possible  without  serious  overlap  of  nearby  argon  lines. 

Initially,  the  horizontal  micrometer  translation  stage  was  adjusted  so  that  the  inlet 
to  the  fiber  optic  cable  was  located  directly  off  the  tip  of  the  cathode  in  the  image  plane. 
From  this  location  on  the  cathode  axis,  the  image  plane  could  be  translated  vertically,  so 
that  the  fiber  inlet  moved  radially  within  the  plasma  plume.  The  vertical  micrometer  was 
adjusted  in  increments  of  0.083  in  (2.108  mm)  between  readings,  which  were  taken  in 
sets  of  3 1 .  One  reading  was  taken  on  the  cathode  axis,  and  15  readings  were  taken  to 
either  side,  up  or  down,  of  the  cathode.  Thus,  a  set  of  3 1  readings  covered  a  "sliver"  of 
the  plasma  jet  2.49  in  (63.246  mm)  in  diameter,  centered  about  the  cathode  axis. 

1.         Data  Format 

As  stated  above,  a  10  A  interval  was  scanned  about  the  location  of  the  Hp 

o  o 

(4861.33  A)  line,  from  4856  to  4866  A  at  a  resolution  of  500  data  points  per  angstrom. 
Each  scan  was  stored  on  a  Macintosh  personal  computer,  using  LabView  software,  as  a 
file  consisting  of  two  arrays  -  wavelength,  recorded  in  angstroms,  and  intensity,  recorded 
as  a  voltage  level.  Each  array  consists  of  5000  double  precision  floating  point  data 
values.  For  bookkeeping  purposes,  each  set  of  scans  was  given  a  group  name,  and 
operating  conditions  such  as  argon/hydrogen  flow  rates,  pressure,  voltage  and  current 
were  also  recorded  for  each  set  of  data.  The  cathode  temperature  profile  was  also 
monitored  during  data  collection. 
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Figure  4.2.  Scan  Location  With  Respect  to  Cathode  Tip. 
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V.  ANALYSIS 

In  order  to  determine  electron  temperatures  and  denisities  from  the  scanned  line 
profiles,  a  computer  program  to  determine  characteristic  line  widths  (both  Doppler  and 
Stark)  was  written  in  Interactive  Data  Language  (DDL).  Sequentially  reading  all  31  files 
in  a  data  set,  the  program  attempts  to  curve  fit  a  Voigt  profile  to  each  line  in  order  to 
estimate  the  parameters  needed  to  extract  the  appropriate  line  widths.  The  appropriate 
characteristics  are  then  determined  either  by  a  simple  formula  in  the  case  of  temperature, 
or  by  interpolation  of  tabulated  theoretical  data  in  the  case  of  electron  densities.  Three 
separate  sets  of  data  were  analyzed,  varying  in  distance  from  the  cathode  tip  and 
operating  parameters,  for  which  an  overview  of  the  analysis  is  given  in  this  chapter.  A 
documented  listing  of  the  DDL  program  is  presented  in  Appendix  A,  and  tabulated  results 
of  the  analysis  are  presented  in  Appendix  B. 
A.         ANALYSIS  PROGRAM 

The  analysis  program  consists  of  eleven  procedures  which  perform  the  tasks  of 
reading  and  preparing  the  data,  determining  the  contributions  of  Doppler  and  Stark  effects 
to  line  width,  calculating  the  electron  temperatures  and  densities,  and  storing  the  results 
in  a  single  data  file  for  a  set  of  scans.  The  majority  of  the  procedures  are  used  in  the 
fitting  of  the  data  and  determination  of  line  widths.  The  program  is  initiated  with  a  single 
call  to  a  master  procedure,  called  LINBROAD,  which  executes  the  above  described 
functions  in  order.  In  turn,  LINBROAD  executes  three  sub-procedures: 

1 .  HTHLE,  which  fits  the  appropriate  Voigt  profile  to  the  scan  and  determines 
the  Stark  and  Doppler  widths  of  the  line  profile. 

2.  NETE,  which  takes  the  line  widths  output  from  HTHLE  and  computes  the 
electron  temperature  and  density  of  the  profile. 

3.  STORIT,  which  writes  the  results  to  a  data  file  with  a  user-specified  name. 
The  first  two  procedures  in  turn  call  one  or  more  sub-procedures  of  their  own  to 

accomplish  their  tasks. 
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1.  FITFILE 

FTTFTLE  is  an  iterative  procedure  which  first  reads  each  data  file  one  at  a  time 
using  a  procedure  called  GETFILE.  Following  the  input  of  a  single  data  file,  FITFILE 
successively  fits  Voigt  profiles  to  the  data,  adjusting  the  guessed  Doppler  width  while 
keeping  track  of  the  Chi-Squared  parameter  of  the  fit,  until  a  minimum  Chi-Squared  is 
reached,  indicating  a  best  fit. 

Once  a  best-fit  Voigt  profile  is  found,  the  HALFWIDTH  procedure  is  called  in 
order  to  determine  the  half-width  of  the  fit.  Following  that,  the  GETPAR  procedure  is 
executed  to  determine  the  fractional  contribution  of  both  Doppler  and  Stark  broadening  to 
the  line  width.  The  PROWIDTH  procedure  is  then  executed,  which  determines  the  full 
half  width  of  the  data.  Lastly,  FITFILE  computes  the  Doppler  and  Stark  widths  of  the 
data  from  the  outputs  of  the  GETPAR  and  PROWIDTH  procedures. 

a.  GETFILE 

GETFILE  first  reads  the  data  file  and  chops  off  the  first  and  last  200  data 
points  in  each  array  (wavelength  and  intensity),  as  the  extreme  wings  of  the  scanned 
profile  contain  significant  portions  of  neighboring  argon  lines  which  would  otherwise 
adversely  affect  the  fitting  of  the  data.  The  procedure  PROFCOR  is  called,  which 
initially  fits  a  Lorentzian  profile  to  the  data  in  order  to  determine  the  location  of  the  line 
center.  PROFCOR  also  determines  an  offset  parameter,  as  the  voltage  (intensity)  levels 
of  the  scan  settle  into  negative  values  out  in  the  line  wings,  which  must  be  corrected  for. 

After  the  intial  fit,  the  profile  is  further  reduced  in  size  to  3000  data  points, 
and  the  profile  spans  3  A  to  either  side  of  the  line  center.  The  profile  is  then  "folded"  in 
half  about  the  line  center,  averaged  and  then  "unfolded"  to  produce  a  smoother,  less 
"noisy"  profile  to  accomodate  the  fitting  of  a  Voigt  profile. 

b.  GETPAR 

Once  a  Voigt  profile  has  been  fitted  to  the  data  and  the  half  width  of  the 
Voigt  profile  is  determined,  GETPAR  is  called  to  determine  the  relative  contributions  of 
Doppler  and  Stark  broadening,  based  upon  the  Doppler  width  and  half  width  of  the  fitted 
Voigt  profile.  This  is  done  by  executing  a  second-order  polynomial  fit  of  the  ratio  of  the 
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Lorentzian  half  width  to  the  full  Voigt  half  width  as  a  function  of  the  ratio  of  Doppler 
half  width  to  the  full  Voigt  half  width.  The  data  to  be  fitted  is  presented  in  Huddlestone 
and  Leonard  (1965),  and  is  reproduced  in  Appendix  C. 
c.  PROWIDTH 

PROWIDTH  determines  the  full  half  width  of  the  observed  profile  by  first 
dividing  the  data  into  bins  of  20  points  each  and  taking  the  average  of  those  bins, 
reducing  the  profile  from  3000  to  150  points,  and  then  finding  the  location  of  the  bin  on 
either  side  of  the  line  center  that  most  closely  represents  the  half  width.  The  right  and  left 
side  half  widths  are  then  added  to  produce  the  full  half  width  of  the  profile  which,  along 
with  the  ratios  determined  in  GETPAR,  is  used  to  compute  Doppler  and  Stark  full  half 
widths. 

2.  NETE 

Given  values  for  the  Doppler  and  Stark  widths  of  the  line  profile,  NETE 
determines  first  the  temperature,  and  then  the  electron  density  of  the  plasma.  The 
temperature  determination  is  easily  accomplished  by  manipulating  equations  (3.4)  and 
(3.15),  giving 

A 1 1/2  2 


T    =  — — K  (5.1) 
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where  AAD1/2  is  the  Doppler  half  width. 

Electron  densities  are  determined  by  interpolation  of  a  table  of  theoretical  Stark 
widths  (in  log  space)  as  a  function  of  temperature,  as  presented  in  Appendix  EI  of  Griem 
(1974).  The  portion  of  this  table  used  is  presented  in  Appendix  C.  The  Stark  widths  for  a 
given  electron  density  as  a  function  of  electron  temperature  are  fitted  with  a  cubic  spline, 
which  is  then  interpolated  for  the  temperature  determined  in  (5.1).  The  spline  is 
performed  in  log  space  with  respect  to  temperature,  resulting  in  the  log  of  the  electron 
density  as  a  linear  function  of  the  log  of  the  Stark  width.  A  linear  interpolation  is  then 
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performed  for  the  log  of  the  observed  Stark  width  of  the  profile.  The  inverse  logarithm  of 
this  interpolated  value  yields  the  electron  density  for  the  profile. 

3.         STORIT 

After  the  electron  densities  and  temperatures  have  been  determined  for  a  data  set, 
STORIT  writes  the  density  and  temperature  vectors,  as  well  as  the  full  half  width, 
Doppler  and  Stark  widths,  and  the  ratios  of  the  Doppler  and  Stark  widths  to  the  full  half 
width  to  a  file  with  a  user-specified  name.  STORIT  also  records  the  distance  of  each 
measurement  from  the  cathode  axis. 

B.  DATA  REPRESENTATION 

Each  set  of  data  is  given  a  name  for  identification,  consisting  of  the  letters  "CTF\ 
for  the  cathode  test  facility  at  JPL,  followed  by  a  number  which  designates,  in  sequence, 
the  running  of  the  CTF  for  that  purpose.  For  example,  the  data  set  CTF144  denotes  the 
144th  running  of  the  facility  for  gathering  data.  The  three  data  sets  analyzed  here  are 
CTF's  144,  145,  and  146. 

For  each  data  set,  operating  conditions  which  may  have  some  bearing  on  the 
analysis  results  are  recorded.  These  are  argon  and  hydrogen  flow,  in  SLPM,  tank 
pressure  in  Torr,  and  operating  voltage  and  current. 

Scans  within  a  data  set  are  numbered  sequentially,  from  1  to  31,  with  the  file 
name  "scanl",  "scan2",  etc.,  "scanl"  being  located  on  the  cathode  axis.  Scans  2  through 
16  make  up  the  15  scans  above  the  cathode,  and  scans  17  through  31  lie  below. 

C.  DATA  ANALYSIS  RESULTS 

Of  the  three  sets  of  data  analyzed,  one,  CTF144,  was  taken  at  a  location  just  off 
the  tip  of  the  cathode,  while  the  remaining  two  were  taken  at  a  distance  of  approximately 
one  half  inch  (1.27  cm)  from  the  cathode  tip.  Tank  pressure,  hydrogen  to  total  volume 
flow  ratio,  voltage  and  current  were  similar  for  CTF  144  and  CTF  145,  with  CTF  146  being 
taken  at  slightly  higher  values.  While  differing  in  numerical  values,  electron  density  and 
temperature  behavior  across  the  plasma  plume  were  similar  for  all  three  data  sets.  Table 
5.1  summarizes  the  operating  conditions  and  extreme  values  for  electron  density  and 
temperature  for  all  three  data  sets. 
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Data  Set 

CTF144 

CTF145 

CTF146 

Distance  from 
Cathode 

Ocm 

1.27  cm 
(0.5  in) 

1.27  cm 
(0.5  in) 

Argon/Hydrogen 
Flow  Rate 

6.08  Ar 
0.8    H2 

6.07  Ar 
0.9    H2 

6.12  Ar 
1.0  H2 

Pressure 

12.25  Torr 

11.8  Torr 

21.4  Torr 

Volatge 

17.6  V 

17.8  V 

22.9  V 

Current 

1400  A 

1400  A 

1400  A 

MaxT 
MinT 

8600  K 
6200  K 

7200  K 
5000  K 

15600  K 
12900  K 

MaxNe 
MinNe 

2.63xl014cm3 
2.30xl014cnV3 

2.41xl014cm-3 
2.14xl014cm'3 

3.95xl014cm3 
3.61xl014crn3 

Table  5.1.  Summary  of  Operating  Conditions  and  Analysis. 

Figures  5.1  and  5.2  below  show  the  radial  electron  density  and  temperature 
profiles  of  all  three  data  sets.  All  results  were  obtained  using  LINBROAD.  Radius  from 
the  cathode  axis  is  represented  on  a  linear  scale  and  electron  density  and  temperature  on  a 
logarithmic  scale.  Apart  from  small  fluctuations  in  Ne  and  Te,  all  three  data  sets  are 
reasonably  well  behaved.  They  suggest  that  the  density  and  temperature  distributions 
across  the  cathode  jet  were  nearly  constant.  No  apparatus  presently  exists  at  the  NAS  A- 
JPL  cathode  test  facility  for  direct  measurement  of  electron  densities  and  temperatures 
within  the  vaccuum  chamber,  and  so  these  results  could  not  be  checked  independently. 
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Figure  5.1.  Radial  Electron  Density  Profiles  Obtined  for  CTF's  144,  145  and  146  Using 
the  LINBROAD  Program.  Tank  pressure,  hydrogen  to  total  volume  flow  rate,  cathode  / 
anode  voltage  and  cathode  current  are  listed  for  each  data  set. 
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X-CTF144:  12.3  Torr,  11.6%  H,  17.6  V,  1400  A 
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Figure  5.2.  Radial  Electron  Temperature  Profiles  Obtained  for  CTF's  144,  145  and  146 
Using  the  LINBROAD  Program.  Tank  pressure,  hydrogen  to  total  volume  flow  rate, 
cathode  /  anode  voltage  and  cathode  current  are  listed  for  each  data  set. 

1.  CTF144 

Data  were  first  collected  for  R  >  0.  Then,  the  fiber  optic  inlet  was  centered  again, 
and  data  was  collected  for  R  <  0.  During  the  collection  of  the  lower  half  of  the  data  set, 
operating  conditions  within  the  vacuum  chamber  were  observed  to  change.  An  increase 
of  hydrogen  flow  from  0.8  to  1.4  SLPM  was  observed,  along  with  a  major  upward  shift  in 
the  cathode  temperature  profile  near  the  tip.  This  shift  could  be  responsible  for  the  slight 
change  in  densities  and  temperatures  shown  in  the  region  of  -1  <  R  <  0  cm  in  Figure  5.1. 
It  appears  as  if  electron  densities  and  temperatures  reached  a  sort  of  "steady  state" 
condition  towards  the  end  of  the  recording  of  the  data  set.  Figure  5.3  below  shows  the 
center  scan  (R  =  0  cm)  from  this  data  set,  before  being  analyzed.  Note  the  two  nearby 
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argon  lines,  one  at  4866A  and  the  other,  weaker  line  beginning  to  show  at  around  4856A. 
It  can  be  seen  that  a  fair  amount  of  noise  is  present  in  the  wings  of  the  486 1.33  A 
hydrogen  line, as  well  as  at  the  line  center.  The  scaling  of  the  intensity  into  both  positive 
and  negative  voltage  values  is  also  evident.  It  was  this  offset  that  necessitated  an  "offset 
constant"  be  included  in  the  curve  fitting  programs. 


FWHM  =  0.85  Angstroms 
T  =  7000  K 
Ne  =  2.35e+14crrvX3 


4856       4857       4858       4859       4860       4861        4862       4863       4864       4865       4866 

Wavelength,  Angstroms 


Figure  5.3.  CTF144  Center  Scan,  Raw  Data. 

2.  CTF145 

On  the  average,  the  CTF145  lines  were  narrower  than  those  in  CTF144,  despite 
similar  operating  voltage  and  current.  Hydrogen  flow  rate  was  slightly  higher  than 
CTF144,  with  a  difference  in  pressure  of -0.5  Torr.  No  shifts  in  operating  conditions 
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were  observed  during  the  data  collection,  and  again  temperature  and  density  behavior  is 
fairly  constant  throughout  the  profile.  Maximum  and  minumum  density  values  differ  by 
12.6%,  and  density  values  are  between  7  and  8%  less  than  those  observed  in  CTF144. 
Temperatures  are  between  16  and  19%  less  than  those  of  CTF144  as  well.  The  decreased 
operating  pressure  during  this  data  collection  may  contribute  to  this  "falling  off  of 
density  and  temperature.  Being  that  the  cathode  jet  is  a  fairly  well  contained  plume  of 
plasma,  both  by  magnetic  and  tank  pressure,  it  is  doubtful  that  the  location  of  this  data  set 
a  mere  one  half  inch  from  the  cathode  tip  would  show  any  major  differences  in  electron 
densities  of  temperatures  than  one  taken  directly  off  the  cathode  tip. 

Figure  5.4  shows  the  center  scan  for  CTF144.  The  realtive  intensity  of  the 
accompanying  argon  lines  is  noticeably  decreased.  The  higher  hydrogen  flow  rate  (12.9% 
of  total  vice  1 1.8%)  may  partially  account  for  this  by  increasing  the  number  of  hydrogen 
emitters  relative  to  argon. 
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Figure  5.4.  CTF145  Center  Scan,  Raw  Data. 
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3.  CTF146 

Scan  30  of  CTF146  could  not  be  completely  loaded  from  disk  due  to  a  disk  error, 
so  the  data  set  was  shrunk  to  27  data  files  by  eliminating  the  last  two  scans  on  either  end 
of  the  data  set.  For  this  set,  tank  pressure  and  voltage  were  significantly  increased,  to 
21.4  Torr  and  22.9  volts  respectively.  The  resulting  line  profiles  are  significantly  broader, 
reflecting  much  higher  temperatures  and  electron  densities  compared  to  the  other  two 
data  sets.  Under  the  increased  pressure  the  electron  density  distribution  is  much  more 
uniform,  with  maximum  and  minimum  values  differing  by  only  9.4%,  with  no  increasing 
or  decreasing  trends  evident  across  either  the  temperature  or  the  pressure  profile.  With 
increased  power,  (VI)  the  hotter  cathode  is  stripping  more  electrons  and  generating  more 
ions,  while  the  increased  pressure  acting  on  the  cathode  jet  forces  the  particles  in  the 
plume  closer  together  into  a  denser,  more  uniform  distribution. 

The  CTF146  center  scan,  shown  in  Figure  5.5  is  very  noticeably  broadened, 
especially  in  the  line  wings  and  the  "transition  region"  between  the  wings  and  line  center. 
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Figure  5.5.  CTF146  Center  Scan,  Raw  Data. 
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D.        ANALYSIS  COMPLICATIONS 

The  scaling  of  the  intensity  values  of  the  line  profiles  presented  problems  in  the 
analyzing  of  the  data,  particularly  in  the  fitting  routines  used  to  model  the  data.  As  is 
evident  in  Figures  5.3,  5.4,  and  5.5,  the  voltages  representing  line  intensities  are  scaled 
such  that  the  line  wings  settle  out  at  negative  levels.  Statistical  distributions,  by  their 
nature,  do  not  accommodate  negative  values,  and  trying  to  curve  fit  a  Voigt  profile  to  one 
of  these  data  sets  is  difficult,  requiring  the  inclusion  of  an  "offset  parameter"  to  estimate 
the  actual  (positive)  background  levels  with  which  the  profile  merges  in  the  line  wings. 
Attempting  to  fit  the  profiles  with  the  LORENTZ  function  written  into  the  program  posed 
little  problem,  outside  of  ambiguities  in  background  intensity  levels  resulting  from  the 
estimation  of  the  offset.  However,  the  VOIGT  function  inherent  to  the  DDL  language  ran 
into  problems  when  attempting  to  fit  profiles  stretching  into  negative  voltages,  and  even 
with  an  offset  parameter  figured  into  the  fit,  the  initial  guesses  for  fit  parameters  had  to  be 
carefully  chosen  to  keep  the  routine  from  producing  meaningless  output. 

Attempting  to  first  offset  the  data  into  positive  voltages  before  attempting  to  fit 
the  data  proved  to  be  useless,  and  no  good  estimate  of  background  levels  could  be  arrived 
at  outside  of  the  offset  parameter  generated  by  curve  fitting  functions.  Though  reasonably 
accurate,  the  fits  generated  for  the  profiles  were  not  as  good  as  expected,  and  some  small 
amount  of  accuracy  was  no  doubt  sacrificed. 


53 


54 


VI.  CONCLUSIONS 

As  stated  above,  no  means  of  directly  measuring  electron  densities  are  in  place  at 
JPL,  and  this  work  can  only  be  considered  a  first  (and  so  far,  only)  estimate  as  to  the 
actual  electron  densities  of  plasmas  generated  in  the  CTF.  The  fact  that  this  technique  is 
known  to  be  the  most  accurate  method  of  non-interfering  analysis  of  plasmas 
(Huddlestone  and  Leonard,  1965),  however,  suggests  that  the  results  produced  by  this 
program  should  provide  a  reasonable  estimate  of  the  actual  values. 

The  results  of  the  three  data  sets  agree  with  what  is  to  be  expected.  Variations  in 
pressure  on  the  plasma  jet  produce  corresponding  variations  in  particle  densities  and 
temperatures;  this  is  particularly  evident  when  comparing  CTF146  to  the  other  two  data 
sets.  Additionally,  the  uniform  radial  electron  density  distribution  agrees  with  Jahn's 
(1968)  statement  that  an  MPD  cathode  jet  has  a  generally  invariant  radial  electron 
distribution. 
A.        FUTURE  TOPICS 

In  addition  to  half  width  measurement,  actual  matching  of  observed  line  profiles 
to  theoretical  Stark  profiles  is  a  proven  method  of  particle  density  determination  with  an 
accuracy  much  greater  than  half  width  matching  alone  (Vidal,  Cooper  and  Smith,  1972). 
A  possible  way  of  doing  this  would  be  by  taking  the  Fourier  transform  of  the  line  profiles 
and  de-convolving  them  with  the  known  Doppler  profile  for  the  plasma.  The  de- 
convolved profile  could  then  be  matched  with  known  line  profiles  to  achieve  an  accurate 
measurement  of  electron  densities.  This  method  requires  accurate  knowledge  of  the 
electron  temperature,  as  well  as  background  intensity  levels  determined  in  a  manner  other 
than  presented  in  this  thesis. 

This  thesis  was  limited  to  three  data  sets  only,  which  were  taken  from  inside  the 
cathode  jet  of  an  MPD  test  facility.  Additional  regions  of  exploration  within  the  plasma 
jets  generated  in  an  MPD  thruster  could  be  near  the  surface  of  the  anode,  within  the  anode 
jet,  in  the  regions  between  the  two  jets,  and  in  the  sheath  region  around  the  cathode.  The 
latter  instance  would  require  a  translating  instrument  with  far  finer  resolution  than  the 
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micrometer  used  in  this  experiment,  being  that  the  sheath  is  on  the  order  of  several 
nanometers  thick. 

Additionally,  a  means  of  direct  measurement  of  electron  densities  and 
temperatures  in  the  vacuum  chamber  should  also  be  explored  as  a  way  of  providing 
independent  verification  for  the  analysis  presented  here. 
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APPENDIX  A.  CODE 

Following  is  the  Interactive  Data  Language  (DDL)  code  developed  to  determine  the 
electron  densities  and  temperatures  for  the  data  sets  taken.  It  consists  of  thirteen  separate 
EDL  procedures,  and  is  run  by  making  a  call  in  IDL  to  a  single  procedure  called 
LINBROAD,  called  as  such: 

LrNBROAD,'<filename>' ,  dx,fwhm,dr,lr,dw,lw,Te,N 
Table  A.  1  defines  the  parameters  as  listed  in  the  above  statement: 


Parameter 

Definition 

<filename> 

Destination  file  for  results. 

dx 

Radial  position  of  scan,  in  mm 

fwhm 

Full  half  width  of  data 

dr 

Ratio  of  Doppler  half  width  to  full 
half  width 

dl 

Ratio  of  Lorentzian  half  width  to  full 
half  width 

dw 

Doppler  width,  in  angstroms 

lw 

Lorentzian  width,  in  angstroms 

Te 

Electron  temperature,  in  K 

N 

Electron  density,  in  cm-3 

Table  A.  1 .  Definition  of  Input/Output  Parameters  for  IDL  Code. 

The  names  of  the  parameters  are  entirely  up  to  the  user,  within  the  restrictions  of  filename 
size  and  restricted  keywords  for  the  IDL  language. 
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;  LINBROAD 

;  Master  procedure  for  analysis  of  broadened  486.133  nm  hydrogen  line 

;  Input  Variables: 

;  ctfno    Name  of  output  file  for  data  storage 

;Output  Variables: 

dx  Position  of  scan  (mm) 

fw  full  half-width  of  line 

dr  Ratio  of  Doppler  half  width  to  observed  half  width  (full) 

lr  Ratio  of  Lorentzian  half  width  to  observed  half  width  (full) 

dw  Doppler  width,  angstroms 

lw  Lorentzian  (Stark)  width,  angstroms 

Te  Electron  temperature,  K 

N  Electron  density,  cmA-3 

pro  linbroad,ctfno,dx,fw,dr,lr,dw,lw,Te,N 

fitfile,lr,dr,fw,dw,lw,dx  ;Call  procedure  FITFILE 

NeTe,fw,lw,dw,Te,N  ;Call  procedure  NETE 

storit,ctfno,dx,fw,dr,lr,dw,lw,Te,N  ;Call  procedure  STORIT 

return  ;Retura 


end 

©fitfile 

@NeTe 
@storit 


;End  procedure 

;Compile  FITFILE 
;Compile  NETE 
;Compile  STORIT 
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;FITFILE 

;  Procedure  sequentially  reads  in  'scanl'  through  'scan31\  preps  and  centers  data,  and 
;  determines  the  Doppler  and  Stark  widths  of  each  observed  line  profile. 

;Output  Variables 

;  lr  Ratio  of  Lorentzian  half  width  to  observed  (full)  half  width 

;  dr  Ratio  of  Doppler  half  width  to  observed  (full)  half  width 

;  fw  Full  half  width  of  observed  profile,  angstroms 

;  dw  Doppler  width  of  observed  profile,  angstroms 

;  lw  Lorentzian  (Stark)  width  of  observed  profile,  angstroms 

;  dx  Position  of  scan  (mm) 


pro  fitfile,lr,dr,fw,dw,lw,dx 

junk  =  check_math(0,0) 

fw  =  dblarr(31) 
dw  =  fw 
lw  =  fw 
lr  =  lw 
dr  =  lr 

dx  =  (findgen(30)  -  15)*0.083*2.54 

for  m  =  0,30  do  begin 

case  1  of 


;Define  output  arrays  fw,  dw,  lw,  lr,  dr 


;Generate  dx 


;Determine  scan  to  read 


(mltl5):  filename  =  strcompress('scan'+string(31-m),/remove_all) 
else  :  filename  =  strcompress('scan'+string(m-14),/remove_all) 


endcase 

getfile,filename,wl,int,par,lfit 
delta  =  sqrt(par(3)) 

lint  =  int-par(O) 

lfit  =  lfit-par(O) 


;Call  procedure  GETFTLE 

;Initial  Doppler  width  estimate  (half  width  of 

;Lorentzian  fit  performed  in  GETFTLE 

; Adjust  scan  and  fit  for  offset  determined  in 

GETFILE 
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a=  fltarr(4) 

a(0)  =  lint(1500) 
a(l)  =  2.75 
a(2)  =  -delta 
a(3)  =  -par(O) 

b  =  a 

w  =  fltarr(n_elements(wl))+l 

chsq  =  999 
chsqp=  1001 

vfit  =  dblarr(3000) 

while  chsq  It  chsqp  do  begin 

chsqp  =  chsq 
pre  =  a(2) 
vpre  =  vfit 


Define     initial-guess     array     for     Voigt 
parameters 

Peak  intensity 

Damping  parameter 

Doppler  width 

Offset  parameter 

Store  initial  guesses 

Define  weight  vector  (no  weighting) 

Set  Chi-Squared  to  999 

Define  previous  Chi-Squared  value  at  1001 

Define  array  to  store  fit  of  Voigt  profile 

Execute  loop  until  Chi-Square  minimized 

Update  chisqp 

Store  current  Doppler  width 

Store  current  fit 


vfit  =  curvefit(wl,int,w,a,function_name  =  Voigt_a',/noderivative)    ;Fit  Voigt  profile 


vint  =  int  +  a(3) 

vfit  =  vfit  +  a(3) 

dop  =  a(2) 

chsq  =  total((vint-vfit)A2/vfit) 

a  =  b 
a(2)  =  dop 

endwhile 

vfit  =  vpre 

halfwidth,wl,vfit,del,ind 
getpar,pre,del,damp,l,d 

lr(m)  =  1 


; Adjust  intensity  by  offset  from  fit 
;Adjust  fit 

;Store  Doppler  width 
;Compute  Chi-Squared  of  fit 

;Reset  a  with  initial  values,  using  Doppler 
;width  from  last  fit 

;Re-execute  fit 

;Best  fit 

;Call  procedure  HALFWIDTH 

;Call  procedure  GETPAR 

;Lorentz  width  fraction 
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dr(m)  =  d 


;Doppler  width  fraction 


prowidth,wl,lint,fwhm 
dw(m)  =  d*fwhm 
lw(m)  =  l*fwhm 
fw(m)  =  fwhm 

print,[dw(m),lw(m),fw(m)] 

endfor 


Call  procedure  PROWIDTH 
Doppler  width 
Lorentzian  (Stark)  width 
Full  half  width 


return 


;Return 


end 

@getfile 
@getdata 
@voigt_a 
@  half width 
@  pro  width 


;End 


;Compile  called  procedures 
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;  GETFILE 


;  Reads  appropriate  scan,  and  fits  an  initial  Lorentzian  profile  to  determine  line  center  and 
;  provide  initial  estimate  parameters  for  fitting  a  Voigt  profile. 

;Input  Variables 

;  filename  Name  of  scan  to  be  read 


;Output  Variables 

;            lam 

Wavelength  vector,  angstroms 

;           int 

Intensity  vector,  volts 

;           par 

Lorentzian  fit  parameters 

;         fit 

Lorentzian  fitted  profile 

pro  getfile,filename,lam,int,par,fit 

wll=dblarr(5000) 

intl=wll 

wlO=wll(0) 

intO=intl(0) 

wl=dblarr(4600) 

val=wl 

st=" 

openr,  1  .filename 

for  j=0,5  do  readf,l,st 

forj=0,4999  do  begin 

readf,l,wlO,intO 
wll(j)=wlO 
intl(j)=intO 
endfor 
close,  1 

wl=wl  1(200:4799) 

val=int  1(200:4799) 

profcor,  wl,val,lam  1  ,int  1  ,coef,fit  1 


;Define  arrays  for  raw  data 

;Set  input  variables  for  reading  data 

;wl  =  wavelength 
;val  =  voltage 
;Blank  string 

;Open  input  file 

;Skip  first  six  lines 

;Read  in  5000  point  wavelength  and  voltage 
vectors 


;Close  input  file 

;Delete  first  and  last  200  points 

;Call  procedure  PROFCOR 
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low  =  where(laml  It  coef(2)-3.0,locount)       ;Set  cutoff  at  +/-  3  angstroms  from  line  center 

high  =  where(laml  gt  coef(2)+3.0,hicount) 

locut  =  locount 

hicut  =  n_elements(laml)-hicount 


lam  =  extrac(laml,locut,hicut-locut) 
int  =  extrac(intl,locut,hicut-locut) 
fit  =  extrac(fitl,locut,hicut-locut) 
par  =  coef 

dl  =  abs(lam(0: 1499)  -  coef(2)) 

dr  =  lam(  1499:2999)  -  coef(2) 

x  =  sort(dl) 

lint  =  int(x(*)) 

rint  =  int(  1499:2999) 

hint  =  (lint+rint)/2 

int(0:1499)  =  hint(x(*)) 
int(1500:2999)  =  hint 

return 
end 

@profcor 


;Extract  working  vectors  from  withing  cutoffs 

;Parameters  from  Lorentzian  fit 
;Split  profile  into  left  and  right  sides 

;"Flip"  left  hand  array 

;Average  left  and  right  hand  arrays 

;"Unflip"  left  hand  array 


;Return 
;End 

;Compile  called  procedures 
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PROFCOR 


Fit  initial  Lorentzian  profile  to  data 


Input  Variables 

lam_0 
int_0 

Wavelength 
Intensity 

Output  Variables 

lam 
int 

coeff 
yfit 

Centered  wavelength 
Returned  intensity 
Fit  coefficients 
Fit 

pro  profcor,lam_0,int_0,lam,int,coeff,yfit 
ain  =  [-1.0,7.0,4861.33,0.21  ;Initial  parameter  guess 


Iorfit,lam_0,int_0,ain,yfit,a 
lam  =  lam_0  -  a(2)+4861.33 

Iorfit,lam,int_0,ain,yfit,a 

int  =  int_0 
coeff  =  a 
return 
end 
@lorfit 


;Call  procedure  LORFIT 

;Set  center  wavelength  to  4861.33  angstroms 

;Fit  around  4861.33  angstroms 


;Return 


;End 


;Compile  called  procedures 
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;LORFIT 


;  Fits  user-defined  Lorentzian  function 


;Input  Variables 


x 

y 

ain 


;Output  Variables 

yfit 

aout 

sigma_a 


X  vector 
Y  vector 
Initial  parameter  guess 


Fitted  Lorentzian 
Returned  parameter  vector 
Standard  deviation  vector 


pro  lorfit,x,y,ain,yfit,aout,sigma_a 
w  =  fltarr(n_elements(x))+l 
aout  =  ain 


;Define  weight  vector  (no  weighting) 


yfit  =  curvefit(x,y,w,aout,sigma_a,function_name  =  'lorentz')        ;Perform  fit 


return 
end 
@  lorentz 


;Return 


;End 


;Compile  fitted  function 
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; LORENTZ 

;Defines  Lorentzian  function 


;Input  Variables 


x 
a 


;Output  Variables 

;         f 

pder 


X  vector 

Parameters  (constants) 


Returned  function 

Array  of  partial  derivatives  wrt  parameters 


pro  lorentz,x,a,f,pder 
f  =  a(0)  +  a(l)  /  ((x  -  a(2))A2  +  a(3)) 
if  n_params()  ge  4  then  begin 

pder  =  fltarr(n_elements(x),4) 

pder(*,0)=  1.0 

pder(*,l)=  1.0/((x-a(2))A2+a(3)) 

pder(*,2)  =  2*a(l)*(x-a(2))/((x-a(2))A2+a(3))A2 

pder(*,3)  =  -a(l)/((x-a(2))A2+a(3))A2 
endif 


;Define  Lorentzian  function 


;Define  partial  derivatives 


end 


;End 
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;  VOIGT_A 

;Define  Voigt  function  for  curve  fitting  based  on  the  DDL  "Voigt"  function 

;Input  Variables 


lam 
a 


;Output  Variables 

;         f 

;  pder 


Wavelength 
Fit  parameters 


Returned  function 

Partial  derivative  vector  (not  used) 


pro  voigt_a,lam,a,f,pder 

u  =  (lam-4861.33)/a(2) 
f  =  a(0)*voigt(a(l),u)-a(3) 

return 

end 


;Dimensionless  wavelength  vector 
;Voigt  function 

;Return 

;End 


67 


;HALFWIDTH 

;  Determines  full  half  width  of  Voigt  profile 

;Input  Variables 


lam 
v 


Wavelength 
Voigt  function 


;Output  Variables 

;  hw 

;  a 


Half  width  of  profile,  in  angstroms 
Index  of  half  maximum  point 


pro  halfwidth,lam,v,hw,a 

val  =  abs((max(v)/2-v)/(max(v)/2)) 
a  =  where(val  eq  min(val)) 

hw  =  abs(lam(1499)  -  lam(a)) 


;Determine  half  max  value 
;Locate  half  max 

;Half  width 


return 


;Return 


end 


;End 
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;GETPAR 

;Determines  damping  ratio,  Lorentzian  and  Doppler  fractional  half  widths  given  doppler 

width 

;and  half  width  of  Voigt  profile. 


;Input  Variables 

;           dw 
;           hw 

Doppler  width 
Half  width 

Output  Variables 

;          d 
;          lr 
;           dr 

Damping  parameter 
Lorentzian  fractional  half  width 
Doppler  fractional  half  width 

pro  getpar,dw,hw,d,lr,dr 

damp  =  fltarr(29) 
lco  =  damp 
dco  =  damp 
dampO  =  damp(O) 
IcoO  =  lco(O) 
dcoO  =  dco(O) 

openr,  1  ,'voigtdat' 

forj  =  0,28  do  begin 

readf,  1  ,damp0,lco0,dco0 
damp(j)  =  dampO 
lco(j)  =  IcoO 
dco(j)  =  dcoO 

endfor 

close,  1 


;Define  input  arrays 
;Define  input  variables 

;Open  input  file 

;Read  Voigt  parameters 


;Close  input  file 
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coef  =  poly_fit(dco,lco,2) 


;Binomial  fit  of  Doppler  and  Lorentz  ratios 


dr  =  abs(dw/hw) 

lr  =  coef(0)+coef(l)*dr+coef(2)*drA2 

d  =  lr/dr 


;Doppler  ratio  of  voigt  profile 

;Lorentz  Ratio 

;  Damping  parameter 


return 
end 


;Return 


;End 
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PROWIDTH 

Determines  full  half  width  of  data.  Profile  width  computed  by  averaging  data  values  over 
a  20 

point  spread,  centered  about  every  data  point  except  the  first  and  last  ten.  Values  are  then 
compared  to  the  average  of  the  20  values  centered  around  the  line  center. 


Input  Variables 


1 


Wavelength 
Intensity 


int 
Output  Variables 

fw  Full  half  width 

pro  prowidth,l,int,fw 

left  =  1(0: 1499) 
right  =  1(1500:2999) 
li  =  int(0:1499) 
ri  =  int(  1500: 2999) 

half  =  total(int(1489:1509))/42 

lcomp  =  dblarr(1480) 
rcomp  =  lcomp 

lval  =  lcomp 
lind  =  lcomp 


rval  =  lcomp 

rind=  lcomp 
is  taken  (right) 

for  a  =  0,1479  do  begin 


;Split  profile  into  left  and  right 


;Average  20  points  around  line  center 
;and  find  half  maximum 

;Comparison    arrays    -    percentage 

difference 

;between  averaged  values  and  half 

max 

;Average  values  of  every  20  points  in  data 

(left) 

;Index  of  point  around  which  average  is  taken 

(left) 

;  Average  values  of  every  20  points  in  data  (right) 
;  Index  of  point  around  which  average 


;Average  20  points  about  every  point  and 

record 

;Index  of  each  point  averaged 
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b  =  a+10 

lval(a)  =  total(li(b-10:b+10))/21 
lcomp(a)  =  abs((lval(a)  -  half)/half) 
lind(a)  =  b 

rval(a)  =  total(ri(b-10:b+10))/21 
rcomp(a)  =  abs((rval(a)  -  half)/half) 
rind(a)  =  b 

endfor 

lx  =  where(lcomp  eq  min(lcomp)) 

lhw  =  abs(left(lind(lx))  -  4861.33) 

rx  =  where(rcomp  eq  min(rcomp)) 

rhw  =  abs(right(rind(rx))  -  4861.33) 

fw  =  rhw+lhw 

return 

end 


;Compute  index  closest  to  half  max 

(left) 

;Compute  half  width  (left) 

;Compute  index  closest  to  half  max 

(right) 

;Compute  half  width  (right) 

;Full  half  width 

;Return 

;End 
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;NETE 

;Computes  electron  density  and  temperature 

;Input  variables 


fw 
lw 
dw 

;Output  Variables 

;  Te 

;  N 


Full  half  width 

Lorentzian  (Stark)  width  (full) 

Doppler  width  (full) 


Temperature  (K) 
Electron  density  (cmA-3) 


pro  NeTe,fw,lw,dw,Te,N 

halfest,fhw,logNe,T,spln 
a  =  n_elements(fw) 
In  =  dblarr(a) 
N  =  ln 

d  =  dw/2 
m=  1.67e-27 
k=1.38e-23 
c  =  3e8 
lamO  =  4861.33 

Te  =  m/(2*alog(2)*k)*(-c*d/lamO)A2 
logTe  =  aloglO(Te) 

hwest  =  fltarr(a,4) 
for  j  =  0,a-l  do  begin 


;Call  procedure  HALFEST 
;Set  counting  variable 
;log(N)  vector 
;N  vector 

;Half  Doppler  width 
;Hydrogen  mass 
;Boltzmann  constant 
;Speed  of  light 
;Center  wavelength 

;Compute  electron  temperature 
;log(Te) 

;2D  array  of  half  width  versus  log(N) 
;to  be  interpolated  from  fhw 

;Spline    interpolation    to    determine 
log(N)  vs 
;log(lw)  for  T 


for  i  =  0,3  do  begin 

hwest(j,i)  =  spl_interp(aloglO(T),fhw(*,i),spln(*,i),logTe(j)) 
endfor 
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lN(j)  =  interpol(logNe,aloglO(hwest(j,*)),aloglO(lw(j)))     ;Linear  interpolation  of  log(N) 

vs 
;log(lw)  to  determine  log(N) 


endfor 
N  =  10A1N 
return 
end 
@halfest 


;Electron  density 

;Return 

;End 

;Compile 
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;HALFEST 


;Read  in  log(N)  vs  alpha  table  and  generate  cubic  spline  for  interpolation  of  electron  density. 
;Output  Variables 


hw 
logN 
T 
spin 


Full  stark  widths  (2D  array) 

Log  (base  10)  of  electron  density  vector 

Electron  temperature  vector 

Spline  coefficients  (2D  array) 


pro  halfest,hw,logN,T,spln 

openr,l,'alpha' 
logN  =  fltarr(4) 
T  =  logN 
hw  =  fltarr(4,4) 
spin  =  hw 

b  =  2.0/3.0 

for  i  =  0,3  do  begin 

for  j  =  0,3  do  begin 

readf,l,T0,N,a0 
T(i)  =  T0*1000 
ifT0eq5thenlogN(j)  =  N 
hw(i,j)  =  a0*2.5e-9*(10AN)Ab 

endfor  ;j 

endfor    ;i 

for  i  =  0,3  do  begin 

spln(*,i)  =  spl_init(aloglO(T),hw(*,i)) 

endfor 

close,  1 

return 

end 


;Open  input  file 

;Define  vectors  and  arrays 


;Electron  density  exponent 

;Fill  T  vector  and  half  width  matrix 


;Generate  splines 


;Close  input  file 

;Return 

;End 
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;STORIT 

;Writes  results  of  analysis  to  file 


;Input  Variables 

;           filename 

File  name  to  write  data  to 

;           dx 

Scan  position,  mm 

;           fw 

Full  half  width 

;           dr 

Doppler  fractional  width 

;          lr 

Lorentzian  fractional  width 

;           dw 

Doppler  width 

;           lw 

Lorentzian  (Stark)  width 

;           Te 

Electron  temperature,  K 

;         n 

Electron  density,  cmA-3 

;Open  write  file 


pro  storit,filename,dx,fw,dr,lr,dw,lw,Te,N 
openw,  1  ,filename 
for  j  =  0,30  do  begin 

printf,l,dx(j),fw(j),dra)4ra)4w(j),lwa),Tea),Na)     ;Write  data  to  file 
endfor 


close,  1 

return 

end 


;Close  write  file 

;Return 

;End 
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APPENDIX  B.  TABULATED  RESULTS 

On  the  following  pages  are,  in  tabulated  format,  some  of  the  results  from  CTF's 
144,  145,  and  146.  Included  are  radial  position,  full  half  width,  Doppler  and  Stark  half 
widths,  electron  temperature,  and  both  the  log  and  decimal  values  for  the  electron  density. 
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Half  Widths,  Angstroms 

R,  cm 

Full 

Doppler 

Stark 

T,  K 

logNe 

Ne,  xlOA14 
cmA-3 

-3.1623 

0.906 

0.3191 

0.7827 

8463 

14.403 

2.529298 

-2.95148 

0.89 

0.296 

0.7816 

7283 

14.404 

2.5351286 

-2.74066 

0.894 

0.2914 

0.7894 

7054 

14.411 

2.5763212 

-2.52984 

0.914 

0.312 

0.797 

8091 

14.415 

2.6001596 

-2.31902 

0.898 

0.3012 

0.7869 

7540 

14.408 

2.5585859 

-2.1082 

0.906 

0.3068 

0.7918 

7823 

14.411 

2.5763212 

-1.89738 

0.902 

0.3157 

0.7808 

8283 

14.402 

2.5234808 

-1.68656 

0.898 

0.3069 

0.7828 

7826 

14.404 

2.5351286 

-1.47574 

0.902 

0.3043 

0.7891 

7694 

14.41 

2.5703958 

-1.26492 

0.906 

0.3217 

0.7808 

8602 

14.401 

2.5176769 

-1.0541 

0.914 

0.3042 

0.8026 

7691 

14.42 

2.630268 

-0.84328 

0.898 

0.3097 

0.7807 

7971 

14.402 

2.5234808 

-0.63246 

0.89 

0.3035 

0.7763 

7653 

14.399 

2.5061093 

-0.42164 

0.85 

0.2787 

0.7494 

6457 

14.379 

2.3933158 

-0.21082 

0.858 

0.2808 

0.7568 

6554 

14.386 

2.432204 

0 

0.85 

0.2894 

0.7417 

6958 

14.372 

2.3550493 

0.21082 

0.842 

0.2848 

0.7361 

6742 

14.367 

2.3280913 

0.42164 

0.85 

0.2965 

0.7365 

7305 

14.366 

2.3227368 

0.63246 

0.846 

0.2974 

0.7314 

7348 

14.362 

2.3014418 

0.84328 

0.842 

0.2852 

0.7358 

6758 

14.367 

2.3280913 

1.0541 

0.842 

0.282 

0.7381 

6609 

14.369 

2.3388372 

1.26492 

0.846 

0.2886 

0.7378 

6922 

14.368 

2.3334581 

1.47574 

0.85 

0.2946 

0.7379 

7214 

14.368 

2.3334581 

1.68656 

0.85 

0.2917 

0.7401 

7070 

14.37 

2.3442288 

1.89738 

0.846 

0.2967 

0.7318 

7317 

14.362 

2.3014418 

2.1082 

0.846 

0.2849 

0.7405 

6744 

14.371 

2.3496328 

2.31902 

0.83 

0.2738 

0.7306 

6229 

14.363 

2.3067472 

2.52984 

0.846 

0.2779 

0.7455 

6419 

14.376 

2.3768403 

2.74066 

0.834 

0.2799 

0.7307 

6511 

14.363 

2.3067472 

2.95148 

0.822 

0.2777 

0.7189 

6409 

14.353 

2.2542392 

3.1623 

0.83 

0.2761 

0.7289 

6337 

14.362 

2.3014418 

Table  B.l.  CTF  144  Results 
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Half  Widths,  Angstroms 

R,  cm 

Full 

Doppler 

Stark 

T,  K 

logNe 

Ne,  xlOA14 
cmA-3 

-3.1623 

0.798 

0.2628 

0.7027 

5740 

14.3405 

2.1902818 

-2.95148 

0.79 

0.2659 

0.6916 

5876 

14.33 

2.1379621 

-2.74066 

0.794 

0.2732 

0.6908 

6201 

14.3282 

2.1291193 

-2.52984 

0.802 

0.2666 

0.7047 

5907 

14.3417 

2.1963422 

-2.31902 

0.814 

0.2666 

0.7178 

5908 

14.3536 

2.2573557 

-2.1082 

0.81 

0.2874 

0.6982 

6865 

14.3333 

2.1542693 

-1.89738 

0.826 

0.2679 

0.7303 

5963 

14.364 

2.3120648 

-1.68656 

0.834 

0.2847 

0.7272 

6735 

14.36 

2.2908677 

-1.47574 

0.834 

0.2813 

0.7297 

6575 

14.3622 

2.3025019 

-1.26492 

0.846 

0.2899 

0.7369 

6985 

14.367 

2.3280913 

-1.0541 

0.834 

0.2864 

0.7259 

6818 

14.3583 

2.2819178 

-0.84328 

0.826 

0.2674 

0.7306 

5942 

14.3648 

2.3163277 

-0.63246 

0.83 

0.2836 

0.7235 

6686 

14.3565 

2.2724796 

-0.42164 

0.822 

0.2832 

0.7149 

6664 

14.3489 

2.233058 

-0.21082 

0.822 

0.2797 

0.7174 

6500 

14.3516 

2.2469841 

0 

0.81 

0.2742 

0.7079 

6250 

14.3438 

2.2069881 

0.21082 

0.822 

0.2834 

0.7147 

6674 

14.3487 

2.2320299 

0.42164 

0.822 

0.2643 

0.7284 

5803 

14.3632 

2.3078097 

0.63246 

0.834 

0.2869 

0.7256 

6842 

14.3579 

2.2798171 

0.84328 

0.83 

0.2975 

0.7265 

6494 

14.3596 

2.2887587 

1.0541 

0.818 

0.2787 

0.7136 

6457 

14.3483 

2.229975 

1.26492 

0.842 

0.2891 

0.7329 

6946 

14.3641 

2.3125972 

1.47574 

0.838 

0.2817 

0.7339 

6595 

14.3658 

2.3216674 

1.68656 

0.834 

0.2888 

0.7242 

6931 

14.3565 

2.2724796 

1.89738 

0.83 

0.2738 

0.7306 

6230 

14.3639 

2.3115325 

2.1082 

0.826 

0.2445 

0.7461 

4966 

14.3814 

2.4065783 

2.31902 

0.826 

0.283 

0.7195 

6655 

14.353 

2.2542392 

2.52984 

0.826 

0.2778 

0.7233 

6413 

14.357 

2.2750974 

2.74066 

0.834 

0.2934 

0.7208 

7153 

14.353 

2.2542392 

2.95148 

0.822 

0.2931 

0.7074 

7140 

14.341 

2.1928049 

3.1623 

0.826 

0.285 

0.718 

6752 

14.3514 

2.2459496 

Table  B.2.  CTF  145  Results. 
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Half  Widths,  Angstroms 

R,  cm 

Full 

Doppler 

Stark 

T,  K 

logNe 

Ne,  xlOA14 
cmA-3 

-2.74066 

1.194 

0.4061 

1.042 

13707 

14.5794 

3.7966451 

-2.52984 

1.206 

0.4162 

1.048 

14393 

14.5827 

3.8256039 

-2.31902 

1.234 

0.4269 

1.072 

15145 

14.5963 

3.9472988 

-2.1082 

1.174 

0.4023 

1.022 

13457 

14.5675 

3.6940264 

-1.89738 

1.206 

0.4188 

1.046 

14577 

14.5814 

3.8141696 

-1.68656 

1.154 

0.3954 

1.005 

12995 

14.557 

3.6057864 

-1.47574 

1.194 

0.4069 

1.042 

13762 

14.579 

3.7931498 

-1.26492 

1.186 

0.4107 

1.03 

14023 

14.5716 

3.7290654 

-1.0541 

1.214 

0.4137 

1.059 

14224 

14.5892 

3.8832916 

-0.84328 

1.21 

0.4078 

1.059 

13820 

14.5894 

3.8850803 

-0.63246 

1.19 

0.4066 

1.037 

13737 

14.5765 

3.7713775 

-0.42164 

1.214 

0.425 

1.051 

15013 

14.5838 

3.8353058 

-0.21082 

1.202 

0.433 

1.031 

15578 

14.5716 

3.7290654 

0 

1.186 

0.3981 

1.039 

13174 

14.5778 

3.7826835 

0.21082 

1.21 

0.4109 

1.057 

14033 

14.5879 

3.8716849 

0.42164 

1.198 

0.4035 

1.048 

13537 

14.5833 

3.8308928 

0.63246 

1.206 

0.4133 

1.05 

14199 

14.584 

3.8370725 

0.84328 

1.194 

0.4073 

1.041 

13788 

14.5788 

3.7914034 

1.0541 

1.182 

0.4133 

1.023 

14199 

14.5676 

3.6948771 

1.26492 

1.182 

0.4025 

1.031 

13462 

14.573 

3.7411059 

1.47574 

1.174 

0.3941 

1.028 

12909 

14.5716 

3.7290654 

1.68656 

1.186 

0.4029 

1.036 

13491 

14.5755 

3.7627035 

1.89738 

1.19 

0.416 

1.03 

14379 

14.5718 

3.7307831 

2.1082 

1.194 

0.416 

1.035 

14373 

14.5746 

3.754914 

2.31902 

1.198 

0.4062 

1.046 

13713 

14.5821 

3.8203223 

2.52984 

1.206 

0.4141 

1.05 

14258 

14.5836 

3.83354 

2.74066 

1.198 

0.4083 

1.045 

13852 

14.5811 

3.8115358 

Table  B.3.  CTF  146  Results 
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APPENDIX  C.  DATA  TABLES 

Two  tables  are  presented  here  which  were  used  in  the  analysis  program  listed  in 
Appendix  A.  Table  C.l,  from  Huddlestone  and  Leonard  (1965)  is  a  tabulation  of  Voigt 
profile  characteristics,  where  XD  is  the  Doppler  half  width,  XL  is  the  Lorentzian  half  width, 
and  Xm  is  the  observed  half  width  of  the  Voigt  profile.  The  first  column,  XL/XD, 
represents  the  damping  parameter  of  the  profile.  The  second  and  third  columns  display, 
respectively,  the  ratios  of  Lorentzian  and  Doppler  half  widths  to  the  profile  width.  Using 
columns  two  and  three,  a  second  order  polynomial  fit  of  XL  vs  XD  was  generated,  from 
which  a  value  of  XL  could  be  determined,  given  XD  and  Xm . 

Table  C.2  from  Griem  (1974)  lists  computed  a1/2  relative  to  electron  densities  and 
temperatures.  These  data  are  interpolated  using  cubic  splines  in  order  to  determine  the 
half  widths  corresponding  to  each  electron  density  for  the  value  of  Te  computed  from  the 
Doppler  half  width  of  the  profile.  Once  the  half  widths  are  determined,  a  linear 
interpolation  (on  a  log-log  scale)  of  Ne  vs  half  width  is  made  for  the  Lorentzian  half  width 
of  the  profile,  and  an  electron  density  is  determined. 
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A, /A, 

A, /A, 

A,/A,. 

oo 

1 

0 

12.01 

0.993 

0.083 

6.01 

0.972 

0.162 

4.00 

0.941 

0.235 

3.00 

0.904 

0.301 

2.70 

0.886 

0.327 

2.40 

0.863 

0.359 

1.80 

0.794 

0.441 

1.50 

0.742 

0.494 

1.20 

0.672 

0.559 

1.14 

0.655 

0.574 

1.08 

0.637 

0.589 

1.02 

0.618 

0.605 

0.96 

0.597 

0.622 

0.90 

0.575 

0.639 

0.84 

0.552 

0.656 

0.78 

0.527 

0.675 

0.72 

0.5 

0.694 

0.66 

0.472 

0.715 

0.60 

0.442 

0.736 

0.54 

0.41 

0.758 

0.48 

0.375 

0.78 

0.42 

0.338 

0.804 

0.36 

0.299 

0.829 

0.30 

0.257 

0.855 

0.24 

0.212 

0.882 

0.18 

0.164 

0.91 

0.12 

0.113 

0.939 

0.05 

0.05 

0.984 

0 

0 

1 

Table  C.l.  Half  widths  of  Voigt  profiles  (From  Huddlestone  and  Leonard,  1965) 
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Jo^Ne^ 

14| 

151           ie|           17I 

Electron  Density,  cm*- 3 

Full  half  width,  alpha  (unitless) 

5000 

7.62e-02 

7.87e-02 

8.08e-02 

7.65e-02 

10000 

7.74e-02 

8.03e-02 

8.40e-02 

8.51  e-02 

20000 

7.77e-02 

8.15e-02 

8.60e-02 

9.02e-02 

40000 

7.90e-02 

8.14e-02 

8.61  e-02 

9.27e-02 

Full  half  width,  angstroms 

5000 

0.4104198 

1 .9675 

9.3760094 

41.203563 

10000 

0.4168831 

2.0075 

9.7473366 

45.835598 

20000 

0.4184989 

2.0375 

9.979416 

48.582502 

40000 

0.4255009 

2.035 

9.99102 

49.929024 

Table  C.2.  Half  width  of  theoretical  Stark  profiles  vs  Ne,  Te.  (From  Griem,  1974). 
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